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STATEMENT  OF  MISSION 


The  mission  of  the  Spectrum  Management  Staff  is  to  assist  the  Department  of 
State,  Office  of  Telecommunications  Policy,  and  the  Federal  Communications 
Commission  in  assuring  the  FAA's  and  the  nation's  aviation  interests  .with 
sufficient  protected  electromagnetic  telecommunications  resources  throughout 
the  world  to  provide  for  the  safe  conduct  of  aeronautical  flight  by  fostering 
effective  and  efficient  use  of  a natural  resource  — the  electromagnetic  radio 
frequency  spectrum. 

This  objective  is  achieved  through  the  following  services: 

• Planning  and  defending  the  acquisition  and  retention  of 
sufficient  radio  frequency  spectrum  to  support  the  aero- 
nautical interests  of  the  nation,  at  home  and  abroad, 
and  spectrum  standardization  for  the  world's  aviation 
community. 

• Providing  research,  analysis,  engineering,  and  evalua- 
tion in  the  development  of  spectrum  related  policy, 
planning,  standards,  criteria,  measurement  equipment, 
and  measurement  techniques. 

• Conducting  electromagnetic  compatibility  analyses  to 
determine  intra/inter-system  viability  and  design  param- 
eters, to  assure  certification  of  adequate  spectrum  to 
support  system  operational  use  and  projected  growth  pat- 
terns, to  defend  aeronautical  services  spectrum  from  en- 
croachment by  others,  and  to  provide  for  the  efficient 
use  of  the  aeronautical  spectrum. 

• Developing  automated  frequency  selection  computer  pro- 
grams/routines to  provide  frequency  planning,  frequency 
assignment,  and  spectrum  analysis  capabilities  in  the 
spectrum  supporting  the  National  Airspace  System. 

• Providing  spectrum  management  consultation,  assistance, 
and  guidance  to  all  aviation  interests,  users,  and  pro- 
viders of  equipment  and  services,  both  national  and  in- 
ternational . 
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A SIMULATION  MODEL  FOR  ANALYZING  PERFORMANCE  OF  SOME 
SIMPLE  COMMUNICATION  SYSTEMS  ON  A DIGITAL  COMPUTER 

Hiroshi  Akima  and  Arthur  D.  Spaulding* 


A d ig i ta 1 -computer  simulation  model  for  analyzing  performances 
of  some  simpie  communication  systems  is  being  developed  for  the  pur- 
pose of  studying  the  effects  of  interfering  signals,  noise,  and/or 
distortions  on  various  communication  systems.  Development  of  the 
model  is  largely  based  on  the  analogy  with  the  laboratory  tests  of 
communication  system  performances.  The  model  consists  of  computer 
subprograms,  each  of  which  either  simulates  a basic  component  of 
communication  systems  or  calculates  characteristics  of  a component. 
An  introductory  explanation  of  computer  simulation,  guidelines  for 
developing  a simulation  model,  an  outline  of  the  model,  and  some 
example  programs  are  presented.  Complete  Fortran  listings  and  a 
user  manual  are  also  presented. 

Key  Words  and  Phrases:  Communication  system,  communication  system 

performance,  computer  simulation,  digital  computer,  distortion, 
interference,  noise. 


1.  INTRODUCTION 

For  efficient  use  of  radio -frequency  spectrum  it  is  important  to  know 
how  communication  systems  perform  when  desired  signals  are  received  with  in- 
terfering signals  and/or  noise  or  suffer  from  distortions  due  to  excessive 
band  limiting  or  othe;  causes.  Analyzing  system  performances  in  such  cases 
is  not  always  easy,  partly  because  demodulators  in  most  systems  are  nonlin- 
ear circuits.  One  way  of  supplementing  our  knowledge  of  system  performances 
is  through  simulation  of  communication  systems  on  a digital  computer.  The 
Institute  for  Telecommunication  Sciences  of  the  Office  of  Telecommunications 
(OT/ITS)  is  developing  a computer  simulation  model  for  this  purpose.  This 
report  describes  the  efforts  for  developing  the  model  and  the  characteristics 
of  the  resulting  model  as  well  as  necessary  information  for  the  use  of  the 
model . 

The  idea  of  computer  simulation  of  communication  systems  is  not  new.  As 
an  example,  we  can  mention  a computer  simulation  model  described  in  CCIR  Re- 
port 520  (CCIR,  1975).  This  CCIR  Report,  however,  does  not  include  program 
listings  of  the  model.  To  the  authors'  knowledge,  no  program  listings  of  sim- 
ulation models  have  been  published  in  the  open  literature. 

According  to  Webster's  dictionary,  the  verb  "simulate"  means  "to  have 
the  external  characteristics  of"  or  "to  act  like".  Therefore,  computer 


*The  authors  are  with  the  Institute  for  Telecommunication  Sciences,  Office  of 
Telecommunications,  U.S.  Department  of  Commerce,  Boulder,  Colorado  80302. 
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simulatioii  cf  a communication  system  is  to  let  the  computer  have  the  external 
characteristics  of  or  act  like  a communication  system.  The  input  data  to  be 
given  to  the  computer  can  be  in  several  forms  — time-sampled  voltage  values 
of  the  input  signal,  frequency-spectrum  components  of  the  input  signal,  sys- 
tem parameters  for  the  modulator,  etc.  In  any  case,  one  must  have  instruc- 
ted the  computer,  through  programming  in  advance,  how  to  act  like  the  com- 
munication system  on  given  data.  The  output  data  from  the  computer  can  also 
take  several  forms  — time-sampled  voltage  values  of  the  output  signal,  fre- 
quencv-spectrum  components,  signal -to-interference  ratio  (SIR)  or  signal-to- 
noise  ratio  (SNR)  values,  etc.  The  simulation  model  is  a major  part  of  the 
instruction  to  be  given  to  the  computer,  and  a well-designed  model  greatly 
reduces  the  user's  difficulty  in  performing  computer  simulation.  (The  reason 
that  the  simulation  model  cannot  be  the  complete  instruction  will  be  explained 
later.) 

In  this  report,  an  introductory  explanation  of  computer  simulation  of 
communication  system  performance  is  given  in  section  2;  the  parallel  rela- 
tionship between  the  computer  simulation  and  the  laboratory  tests  is  empha- 
sized. In  section  3,  guidelines  for  developing  the  model  are  derived.  In 
section  4,  a general  description  of  the  model  developed  in  this  study  is 
given.  Section  5 describes  how  to  use  the  model;  some  simple  examples  are 
included  to  illustrate  its  use.  Section  6 summarizes  all  preceding  sections 
and  gives  some  suggestions  for  further  development,  Fortran  listings  of  all 
subprograms  included  in  the  model  are  given  in  Appendix  A,  and  a user  manual 
is  given  in  Appendix  B. 


2.  SIMULATION  OF  SYSTEM  PERFORMANCE  ON  A DIGITAL  COMPUTER 
2.1.  Need  for  Computer  Simulation 

Typically,  there  are  three  methods  of  studying  performances  of  communi- 
cation systems  including  the  performance  of  demodulators;  i.e.,  theoretical 
analysis,  laboratory  test,  and  computer  simulation.  Each  has  its  advantages. 
Theoretical  analysis  is  convincing  if  the  analysis  is  done  without  approxi- 
mation errors.  There  are,  however,  many  cases  in  which  theoretical  analyses 
do  not  look  feasible  even  if  not  altogether  impossible.  A laboratory  test  is 
also  good  if  the  test  is  done  in  a well -control led  situation.  In  many  cases, 
however,  laboratory  tests  lack  flexibility;  changing  a system  parameter  some- 
times necessitates  rebuilding  of  all  of  the  equipment  in  the  test  arrangement. 
Laboratory  tests  can  also  be  quite  expensive  and  time  consuming,  especially 
tests  for  designing  new  systems.  Computer  simulation  can  mitigate  the  disad- 
vantages of  the  other  two  methods;  it  is  almost  always  possible  (at  least  the- 
oretically if  sufficient  time  is  allowed  for  computation)  and  it  is  flexible. 

As  a typical  example  we  will  take  a frequency-modulation  (FM)  demodula- 
tor and  first  consider  the  case  where  an  FM  desired  signal  is  received  with 
white  Gaussian  noise.  The  output  signal-to-noise  ratio  (SNR)  has  been  pre- 
viously calculated  as  a function  of  input  SNR  and  the  modulation  index  (i.e., 
the  ratio  of  peak  frequency  deviation  to  maximum  frequency  of  the  baseband 
signal)  (Stumpers,  1948;  Akima,  1963).  Second,  we  consider  the  case  where 
the  desired  signal  is  received  with  non-Gaussian  noise.  We  can  even  consider 
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the  third  case  where  the  desired  signal  is  received  with  an  interfering  sig- 
nal and  noise  (either  Gaussian  or  non-Gaussian) . Very  little  is  known  about 
the  output  SNR  for  these  cases.  Solving  even  these  simple  problems  (e.g., 
calculating  the  output  SNR)  by  theoretical  analysis  is  too  complicated. 

We  may  ask:  "Can  we  solve  the  same  problems  with  laboratory  tests?" 

The  answer  is  "Yes,"  of  course,  but  we  need  some  further  considerations.  As 
mentioned  earlier,  the  test  set-ups  are  not  simple  even  for  these  simple  prob- 
lems. If  we  design  all  the  components  such  as  the  modulator,  noise  generator, 
bandpass  filter,  and  demodulator  to  be  flexible  enough  for  possible  changes 
in  system  parameters,  the  cost  will  be  very  high.  If  we  disregard  the  sig- 
nificance of  flexibility  in  initial  test  set-ups,  changing  the  modulation  in- 
dex, for  example,  may  necessitate  rebuilding  of  the  modulator,  bandpass  fil- 
ter, and  the  demodulator. 

This  is  a good  example  of  one  way  in  which  computer  simulation  can  do  a 
better  job.  As  will  be  explained  later,  we  can  write  several  computer  sub- 
programs that  simulate  system  components  such  as  the  FM  modulator,  bandpass 
filter,  ,-M  demodulator,  etc.  If  properly. written,  a computer  subprogram  is 
flexible  in  nature.  A subprogram  that  simulates  an  FM  modulator,  for  example, 
works  for  any  legitimate  values  of  system  parameters  such  as  the  baseband 
bandwidth,  center  frequency,  modulation  index,  etc.  Writing  a flexible  sub- 
program does  not  usually  require  extra  effort  or  extra  cost,  both  of  which 
are  required  for  building  a flexible  hardware  system  component.  Running  lab- 
oratory tests  with  various  sets  of  parameter  values  is  equivalent  to  calling 
the  subprograms  with  various  sets  of  parameter  values.  Computer  simulation 
is  very  useful  when  system  performance  is  studied  for  system  parameters  of 
wide  variability. 

In  the  laboratory  tests,  care  must  be  taken  so  that  test  conditions  will 
not  change  inadvertently.  For  example,  stabilizing  the  power  supply  voltage 
is  required  in  many  cases.  Room  temperature  can  be  another  unstable  condi- 
tion. All  parameter  values  except  those  which  one  wants  to  change  intention- 
ally must  be  kept  constant.  Computer  simulation  takes  care  of  this  require- 
ment automatically.  In  computer  simulation,  only  those  parameters  one  wants, 
to  change  are  changed,  and  other  parameters  remain  constant. 

When  one  studies  the  effects  of  an  interfering  signal  of  a specified 
type,  for  example,  all  test  equipment  must  be  free  from  other  unintentional 
interferences  and/or  noise  (either  electromagnetic,  acoustic,  or  mechanical). 
In  other  words,  intensities  of  unintentional  interference  and  noise  must  be 
negligibly  small  compared  with  those  of  the  desired  signal  and  of  the  inter- 
fering signal  in  question.  Computer  simulation  can  take  care  of  this,  too, 
very  easily.  The  main  source  of  "noise"  in  computer  simulation  is  the  one 
due  to  errors  in  binary  approximation  of  real  numbers,  but  it  is  usually  very 
small.  If  one  simulates,  for  example,  a back-to-back  test  of  an  FM  system  in 
whicn  an  FM  demodulator  is  connected  to  an  FM  modulator  directly,  an  output 
SNR  of  200  dB  or  higher  is  obtained  on  a CDC-6000  series  computer. 


Computer  simulation  was  first  considered  as  a supplement  to  laboratory 
tests  in  the  beginning  of  this  section,  but  the  last  two  paragraphs  indicate 
that  computer  simulation  can  do  more  than  laboratory  tests  in  certain  aspects 
One  can  simulate  an  idealized  laboratory  test  very  easily. 


2.2.  Computer  Simulation  of  System  Performance 


Computer  simulation  of  a communication  system  allows  the  computer  to  act 
like  a communication  system  or  like  a part  of  the  system.  Assuming  a typical 
example,  we  will  first  describe  what  will  take  place  in  the  computer  simula- 
tion. 

Suppose  one  is  studying  the  effects  of  an  AM  interfering  signal  on  an  FM 
desired  signal  in  the  presence  of  white  Gaussian  noise.  The  computer  produces 
data  (a  sequence  of  numbers)  that  represent  the  FM  desired  signal  with  speci- 
fied system  parameters  such  as  the  frequency  of  the  modulating  signal.,  the 
modulation  index,  and  the  center  frequency  of  the  modulated  signal.  The  data 
can  be  a sequence  of  time-sampled  voltage  values  or  a frequency  spectrum  (i.e., 
Fourier  transform  of  the  time  sequence).  (The  choice  as  to  whether  a time 
sequence  or  its  spectrum  is  used  is  a matter  of  technicality;  it  is  discussed 
in  detail  later.)  Similarly,  the  computer  generates  data  for  th;j  AM  interfer- 
ing signal  and  data  for  the  noise.  Next,  the  computer  superposes  these  three 
sets  of  data  with  relative  ratios  that  correspond  to  a specified  signal-to- 
interference  ratio  (SIR)  and  a signal-to-noise  ratio  (SNR);  the  resulting  data 
represent  the  input  signal  to  the  receiver.  Next,  the  computer  acts  like  an 
intermediate-frequency  (IF)  bandpass  filter;  it  modifies  the  receiver  input 
signal  and  obtains  the  filter  output  signal.  Then,  the  computer  acts  like  an 
FM  demodulator;  it  produces  a baseband  signal  that  corresponds  to  the  demodu- 
lator output.  If  so  instructed,  the  computer  does  further  baseband  process- 
ing; it  acts  like  a baseband  filter  or  a psophometer.  (A  psophometer  is  an 
apparatus  for  the  objective  measurement  of  baseband  noise>  it  is  essentially 
a bandpass  filter  that  represents  average  frequency  response  of  human  ears 
and  the  frequency  response  of  a typical  electro-acoustic  transducer  in  the 
receiver.  The  psophometrically-weighted  noise  power  is  usually  a good  mea- 
sure of  the  disturbing  effect  of  the  baseband  noise.)  The  computer  can  even 
act  as  a monitoring  device  such  as  an  oscilloscope  or  a spectrum  analyzer;  it 
can  display  a signal  at  any  point  in  the  system  as  a time  sequence  or  a spec- 
trum. 


As  described  in  this  example,  the  computer  acts  like  each  system  compo- 
nent along  the  flow  of  the  signal  in  the  system.  (To  simulate  two  system 
components  or  more  that  operate  in  parallel  such  as  the  FM  and  AM  modulators 
and  the  noise  generator  in  this  example,  the  computer  simulates  one  component 
at  a time  in  an  arbitrary  order.)  The  flow  of  data  in  the  computer  simula- 
tion is  shown  schematically  in  figure  1.  It  is  clear  from  this  figure  (and 
also  from  the  preceding  descriptions)  that  signal  flow  in  a computer  simula- 
tion parallels  the  flow  in  a laboratory  test.  This  parallel  relationship 
between  computer  simulation  and  laboratory  tests  is  very  important  in  devel- 
oping the  model  for  computer  simulation  as  well  as  in  using  the  model  for 
solving  specific  problems. 

In  the  development  of  a computer  program  for  system  simulation,  it  is 
advisable  to  set  apart  the  portions  of  the  program  that  correspond  to  iso- 
lated procedures.  The  more  subprograms  one  can  write,  the  simpler  the  main 
program  will  be.  Since  it  is  almost  impossible  to  please  every  user  in  the 
way  of  writing  the  input  and  output  parts  of  the  program,  it  is  desirable 
that  each  procedure  that  simulates  each  system  component  - each  box  in  fig- 
ure 1 — be  written  as  a subprogram,  leaving  it  up  to  the  user  to  determine 
the  main  program  which  makes  use  of  the  appropriate  subprograms. 
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Figure  1.  An  example  of  computer  simulation  arrangement. 
(A  voice  signal  case.) 


3.  GUIDELINES  FOR  DEVELOPING  A SIMULATION  MODEL 
3.1.  Main  Structure  of  Computer  Simulation  Model 

In  the  preceding  section,  we  have  considered  the  parallel  relationship 
between  a computer  simulation  model  and  laboratory  test  facilities  and  we  have 
suggested  that  a computer  simulation  model  of  communication  system  performance 
will  be  a collection  of  computer  subprograms,  each  simulating  a system  compon- 
ent. In  this  section,  v<e  proceed  with  the  parallel  relationship  and  determine 
the  main  structure  of  the  computer  simulation  model. 

Suppose  a coinmunication  engineer  is  studying  in  the  laboratory  the  per- 
formance of  a communication  system  in  which  the  desired  signal  is  received 
in  the  presence  of  an  interfering  signal  and  noise.  The  first  thing  he  does 
when  he  enters  the  laboratory  is  perhaps  to  locate  necessary  system  compon- 
ents such  as  the  modulator,  demodulator,  noise  generator,  etc.,  on  the  shelves 
of  the  laboratory.  If  he  can  locate  necessary  components,  he  brings  them  to 
the  test  bench.  He  must  buy,  borrow,  or  build  the  missing  components,  if  any. 
Next,  he  makes  necessary  cable  connections  among  these  components.  To  do  this 
and  also  to  conduct  the  test  that  follows,  he  usually  needs  a user  manual  for 
each  component.  Before  starting  the  test,  he  must  consider  and  decide  many 
things.  The  ranges  and  increments  of  the  signal-to-interference  ratio  (SIR) 
and  signal-to-noise  ratio  (SNR),  the  mjdulation  indexes,  the  frequency  detun- 
ing of  the  signals  are  only  a few  of  such  items.  Decisions  on  these  items  is 
his  responsibil ity,  not  of  the  laboratory  facilities.  The  facilities  usually 
include  various  items  of  test  equipment  such  as  voltmeters,  power  meters,  VU 
(volume  unit)  meters  for  voice  signals,  CRT  (cathode-ray-tube)  oscilloscopes, 
etc.,  and  the  selection  of  these  pieces  of  equipment  for  a specific  test  is 
also  his  responsibility.  During  the  test,  he  may  have  to  modify  his  test 
plan,  depending  on  the  partial  results  of  the  test.  Finally,  he  tries  to 
present  the  test  results  in  a comprehensive  manner.  The  decision  on  the 
means  of  presentation  (i.e.,  tables,  graphs,  empirical  formulas,  etc.)  or  on 
the  format  of  the  tables  or  graphs  is  also  his  responsibility.  This  example 
suggests  several  guidelines  for  developing  laboratory  test  facilities  for 
communication  system  performances. 
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The  ubuve  example  indicates  clear  distinction  between  the  responsibility 
of  the  laboratory  and  that  of  the  user.  The  laboratory  including  its  mainte- 
nance staff,  if  any,  does  not  conduct  a test  by  itself.  All  details  of  plan- 
ning and  conducting  the  test  and  of  representing  the  test  results  are  the  re- 
sponsibility of  the  user.  A good  laboratory  is  the  one  that  provides  the 
maximum  ease  to  the  user.  Since  the  computer  simulation  model  is  a software 
replica  of  a test  laboratory,  the  model  does  not  perform  simulation  by  itself; 
the  user  must  do  this  by  writing  a computer  program.  A good  computer  simula- 
tion model  is  the  one  that  keeps  the  difficulties  to  the  user  to  a minimum. 

A good  laboratory  must  have  a variety  of  communication  system  components. 
Their  types  and  operating  ranges  must  cover  the  users'  needs.  Each  component 
must  work  properly.  It  must  be  well  documented.  Standardized  cable  connec- 
tors and  power  supplies  are  very  desirable  for  easing  the  test-arrangement 
efforts.  These  are  some  of  the  guidelines  for  developing  good  laboratory 
test  facilities. 

Translation  of  these  guidelines  for  good  laboratory  test  facilities  to 
those  for  a good  computer  simulation  model  is  almost  verbatim.  A good  com- 
puter simulation  model  must  have  a variety  of  computer  subprograms,  each  sim- 
ulating a communication  system  component.  Subprogram  types  and  operating 
ranges  must  cover  the  users'  needs.  (In  general,  however,  satisfying  the  re- 
quirements for  wide  operating  ranges  of  parameters  for  a computer  subprogram 
is  not  as  severe  as  for  a hardware  system  component.)  Each  subprogram  must 
do  exactly  what  it  is  expected  to  do;  it  must  have  been  completely  debugged. 
Each  subprogram  must  be  well  documented  so  that  every  user  can  call  it  with- 
out difficulty.  Procedures  of  calling  these  subprograms  must  be  standardized 
and  consistent  within  the  model.  A model  must  be  written  in  a standardized 
computer  language  such  as  American  National  Standards  Institute  (ANSI)  Stan- 
dard Fortran  (ANSI,  1966)  so  that  the  model  can  be  used  on  most  computers. 


3o2„  Frequency- Domain  vs.  Time-Domain  Approach 

The  analogy  with  the  laboratory  test  has  been  extensively  used  above. 

But  there  are  several  things  to  which  this  analogy  does  not  apply.  The  sub- 
ject of  this  section,  i.e.,  the  discussion  of  frequency-domain  approach  vs. 
time-domain  approach,  is  one  of  them. 

In  laboratory  tests,  a signal  voltage  is  continuous  most  of  the  time. 
Because  of  the  sampling  theorem  (Wozencraft  and  Jacobs,  1967,  ch.  8),  a band- 
limited  signal  can  be  represented  by  a sequence  cf  signal  voltage  values  that 
are  measured  at  sampling  points  equally  spaced  in  time.  (The  sampling  rate 
must  be  equal  to  or  higher  than  the  so-called  Nyquist  rate  that  is  equal  to 
twice  the  maximum  frequency  of  the  signal.)  Representation  of  a signal  by  a 
sequence  of  time-sampled  voltage  values  is  not  the  only  way  of  representing  a 
signal.  It  can  also  be  represented  by  its  Fourier  transform,  or  a sequence 
of  frequency  spectrum  components.  These  two  representations  are  equivalent 
mathematically,  and  transformation  from  one  to  the  other  can  be  done  easily 
on  the  computer. 
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Since  these  two  representations  are  mathematically  equivalent,  the  choice 
of  one  depends  on  users'  convenience.  Problems  associated  with  this  choice 
are  discussed  below.  Let  a filter  (either  a bandpass  or  a lowpass  filter)  be  1 

a first  example.  In  general,  characteristics  of  a filter  are  commonly  ex-  | 

pressed  by  its  frequency  responses;  i.e.,  amplitude  response  (or  attenuation)  I 

and  phase  response  (phase  delay)  versus  frequency.  Filter  characteristics  i 

can  also  be  expressed  by  the  impulse  response,  i.e.,  the  transient  output  j 

waveform  when  a unit  impulse  is  imposed  at  the  filter  input.  Effects  of  a | 

filter  can  be  calculated  either  by  multiplying,  on  a component-by-component  | | 

basis,  the  spectrum  components  of  the  input  signal  by  the  filter  frequency  | 

responses  in  frequency  domain,  or  by  convolving  the  input  signal  with'  the  im- 
pulse response  in  time  domain.  These  two  methods  of  calculation  are,  again, 
equivalent  mathematically,  but  they  are  very  different  in  degree  of  conveni- 
ence of  numerical  calculation.  Frequency  response  of  a filter  substantially 
diminishes  if  the  frequency  deviates  from  the  filter  passband  to  some  extent, 
but  it  takes  a long  time  for  the  impulse  response  to  diminish  substantially. 

In  the  frequency  domain,  a spectrum  component  of  the  output  signal  is  a sim- 
ple product  of  the  respective  spectrum  component  of  the  input  signal  and  the 
frequency  response  of  the  filter,  but  in  the  time  domain,  each  sampled  value 
of  the  output  signal  is  an  integral  (usually  a truncation  of  an  infinite  in- 
tegral) of  the  product  of  the  input  signal  and  the  filter  impulse  response. 

Therefore,  we  can  see  that  the  frequency -domain  representation  is  much  more 

convenient  for  representing  the  characteristics  of  a filter.  j 

Modulators  will  be  considered  next.  To  generate  a double-sideband  (DSB)  i 

amplitude-modulation  (AM)  signal,  either  time  of  frequency  domain  approach  is  ! 

equally  simple  and  straightforward.  For  a single-sideband  (SSB)  signal,  the 
frequency-domain  approach  is  much  simpler  than  the  time-domain  approach;  SSB 
modulation  is  equivalent  to  shifting  the  spectrum  components  of  the  baseband  ' 

signal  (modulating  signal)  to  the  radio  frequency  region.  For  other  types  of  | 

modulations  such  as  phase  modulations  or  frequency  modulations  (FM),  generat-  i 

ing  a signal  as  a function  of  time  is  much  simpler  than  generating  its  spec-  : 

trum  directly. 

As  for  the  demodulators,  the  situation  is  a little  different  except  for 
the  SSB  demodulator,  for  which  the  frequency-domain  approach  is  simpler.  To 
demodulate  a DSB-AM  signal,  the  time-domain  approach  is  simpler.  For  an  FM 
signal,  a method  starting  with  the  spectrum  is  not  only  simpler  but  also  more 
accurate,  as  described  later. 

From  these  observations,  we  can  conclude  that  there  is  no  clear  basis 
for  selecting  one  over  the  other,  insofar  as  the  modulators  and  demodulators  i , 

are  concerned.  In  this  situation,  convenient  representation  of  filter  char- 
acteristics with  the  frequency  response  might  be  a decisive  factor.  Of  | 

course,  it  is  possible  to  use  one  representation  for  some  system  components  j 

and  the  other  representation  for  other  components,  but  we  reject  this  idea.  j 

Such  an  idea  would  confuse  the  prospective  users  of  the  model;  it  is  not  easy  \ 

for  many  users  to  remember  which  representation  is  used  for  the  input  and  out-  j 

put  of  a specific  system  component.  j 
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3.3.  Lower-Frequency  Approach 


Complexity  of  computation  increases  with  the  increase  in  number  of  sam- 
pling points  (or  number  of  frequency  components)  to  be  handled  by  the  model. 
It  is  desirable  to  keep  these  numbers  as  small  as  possible.  To  reduce  the 
numbers,  let  us  review  a simple  rule  that  governs  the  number  of  sampling 
points  (or  the  number  of  frequency  components). 


Suppose  we  wish  to  represent  a voltage  waveform  of  duration  T,  which  is 
a real  function  of  time  t,  and  we  do  this  by  taking  N voltage  values  at  N sam- 
pling points  equally  spaced  in  time.  The  minimum  number  N that  is  necessary 
to  represent  the  exact  waveform  depends  on  the  fineness  of  the  waveform  or  the 
maximum  frequency  component  contained  in  the  waveform.  According  to  the  sam- 
pling theorem,  the  sampled  values  taken  at  the  rate  of  2 fm  can  exactly  rep- 
resent the  waveform  when  the  waveform  does  not  contain  any  frequency  spectrum 
components  higher  than  fm.  Therefore,  if  we  take 


N'  = 2 f„ 


samples  equally  spaced  in  time  for  a waveform  having  a duration  T and  maximum 
frequency  fm,  these  sampled  values  exactly  represent  the  waveform.  The  same 
number  applies  to  the  number  of  frequency  spectrum  components  that  represent 
the  same  waveform.  As  is  well-known,  such  a waveform  can  be  represented  with 
a series  of  sine  and  cosine  functions  having  a fundamental  frequency  1/T  and 
harmonic  frequencies  i/T,  where  i = 2,3, ...  . Since  no  component  exists  that 
is  higher  than  fm,  there  are  N'  = 2 fmT  components  in  the  series;  i.e.,  fmT 
sine  components  and  fmT  cosine  components. 


In  the  preceding  paragraph,  we  have  represented  everything  with  real  num- 
bers. From  now  on,  however,  we  will  use  complex  numbers.  We  represent  the 
waveform  with  a complev  function  of  a real  argument  t;  the  actual  waveform  is 
a real  part  of  this  complex  function.  The  cosine  and  sine  functions  of  an 
identical  angle  are  combined  into  a complex  exponential  function  as 


exp  (jx)  = cos  X + j sin  x. 


Using  complex-number  representations,  we  can  represent  a waveform  of  duration 
T and  of  maximum  frequency  fm  with  either 


N = fmT 


time  sampled  values  or  N frequency  spectrum  components. 


It  is  clear  from  the  above  considerations  that  reducing  N is  equivalent 
to  reducing  the  frequency-time  produce,  fmT.  A simple  way  of  achieving  this 
is  to  lower  fm  as  much  as  possible.  Let  us  consider  an  example.  The  spec- 
trum of  an  AM  signal  modulated  with  a voice  signal  having  the  maximum  fre- 
quency of  3 kHz  extends  from  3 kHz  below  the  carrier  frequency  to  3 kHz  above 
the  carrier  frequency.  If  the  carrier  frequency  is  100  kHz,  fm  is  103  kHz, 

If  we  consider  the  signal  in  an  interval  of  10  ms,  N must  be  at  least  1030. 
But,  if  we  lower  the  carrier  frequency  to  5 kHz,  fm  is  8 kHz,  and  N for  10 
ms  is  80.  In  the  laboratory  tests,  modulating  a carrier  of  5 kHz  with  a 


! 

I 


modulating  signal  of  3 kHz  is  neither  easy  nor  advantageous,  but  is  is  defi- 
nitely good  in  computer  simulation. 

This  lower-frequency  approach  should  not  be  confused  with  the  so-called 
lowpass-equivalent  analysis  that  is  used  for  the  analysis  or  synthesis  of  a 
bandpass  filter.  Our  approach  is  equivalent  to  beating  down  the  frequency 
(similar  to  the  process  used  in  a superheterodyne  receiver).  . 


3o4.  Basic  Parameters  for  Simulation  Model 

One  of  the  most  basic  parameters  in  the  computer  simulation  model  is  the 
number  of  frequency-spectrum  components  (or  the  number  of  sampling  points). 
This  number  dictates  the  necessary  core-memory  size  in  the  computer  and  de- 
termines the  computation  time  required  for  simulation.  Since  this  number  is 
the  product  of  the  maximum  frequency  of  the  signal  and  the  time  duration  of 
the  signal,  the  latter  two  are  also  important  parameters. 

The  frequency  of  the  signal  in  the  computer  simulation  model  does  not 
have  to  be  equal  to  the  frequency  of  the  signal  to  be  simulated.  As  in  the 
laboratory  tests,  one  can  choose  an  arbitrary  frequency  as  the  center  fre- 
quency of  the  signal,  as  long  as  the  relative  relations  among  the  center  fre- 
quencies of  the  signals  and  the  bandpass  filter  are  kept  unchanged  and  the 
minimum  frequencies  of  all  the  signals  are  above  zero.  The  maximum  frequency 
used  in  the  model  must  be  equal  to  or  greater  than  the  maximum  frequencies  of 
all  the  signals  involved. 

Representing  a waveform  of  duration  T with  frequency  spectrum  components 
implies  that  we  are  considering  a waveform  of  infinite  duration  consisting  of 
infinite  repetitions  of  the  waveform  of  duration  T.  Even  if  the  model  uses 
the  time-domain  approach  and  represents  the  input  and  output  voltage  wave- 
forms with  sequences  of  time-sampled  voltage  values,  use  of  frequency  spec- 
tra of  these  waveforms  cannot  be  avoided.  Rigorously  speaking,  therefore, 
computer  simulation  of  the  type  we  are  considering  can  handle  a waveform  only 
when  the  waveform  is  a periodic  function  of  T. 

If  we  take  the  frequency-domain  approach,  the  fundamental  frequency  or 
the  unit  frequency  spacing  of  the  spectrum,  fo,  is  an  important  parameter. 

It  is  the  reciprocal  of  T ana  is  equal  to  fm/N.  If  we  take  the  time-domain 
approach,  the  unit  time  interval  of  sampling,  to,  is  an  important  parameter. 
It  is  the  reciprocal  of  fm  and  is  equal  to  T/N. 

•So  far  we  have  identified  five  parameters;  i.e.,  the  number  of  frequency 
spectrum  components  (or  the  number  of  sampling  points)  N,  the  unit  frequency 
spacing  in  the  spectrum  fo,  the  maximum  frequency  fm,  the  unit  time  interval 
of  sampling  to,  and  the  time  duration  of  the  signal  T.  These  five  parameters 
have  their  own  physical  meanings,  but  they  are  not  independent  of  each  other. 
It  is  very  easy  to  show  that,  if  we  select  a pair  of  independent  parameters, 
the  remaining  three  parameters  can  be  represented  as  simple  functions  of  the 
former  two.  We  will  select  two  parameters  for  the  model  in  the  paragraph 
that  follows. 
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Since  N is  the  basic  parameter  in  computer  programming^  we  will  select  N 
as  a basic  parameter  for  the  model.  For  representing  the  characteristics  of 
modulators  and  demodulators,  any  one  out  of  the  remaining  four  parameters  can 
be  used  as  the  basic  parameter.  For  representing  the  characteristics  of  a 
filter,  however,  neither  to  nor  T has  a good  and  comprehensive  meaning.  As 


for  the  choice  of  one  out  of  fo  and  fm,  fo  is  more  convenient  for  represent- 
ing other  frequencies  such  as  the  center  frequencies  of  signals  or  the  band- 


widths  of  filters.  In  summary,  we  will  select  the  number  of  frequency  spec- 


trum components  N and  the  unit  frequency  spacing  in  the  spectrum  fo  as  basic 


parameters  for  the  computer  simulation  model. 


3.5.  Programming  Language 

In  general,  a Fortran  language  is  accepted  by  most  computer  system  in- 
stallations today.  A Fortran  language  is  one  of  the  basic  computer  languages 
taught  in  beginning  programming  classes.  Learning  Fortran  programming  does 
not  take  a long  time  for  average  scientists  and  engineers.  Perhaps  the  only 
tedious,  but  not  difficult,  part  of  Fortran  programning  is  specifying  the 
proper  format  of  data  in  conjunction  with  the  input  and  output  statements. 

There  are  several  versions  (or  dialects)  of  Fortran  that  particular  com- 
puter systems  can  accept.  Conversion  from  one  version  of  Fortran  to  another 
is  sometimes  difficult.  For  the  purpose  of  minimizing  the  conversion  prob- 
lems to  a specific  version,  ANSI  standard  Fortran  (ANSI,  1966)  has  been  de- 
veloped. Most  Fortran  compilers  can  compile  programs  that  are  written  in 
ANSI  standard  Fortran.  When  revisions  of  programs  written  in  ANSI  standard 
Fortran  are  required  by  a computer  system,  such  revisions  are  usually  minor. 
In  such  a case,  the  user  manual  of  the  computer  system  usually  gives  instruc- 
tions on  how  to  revise  the  programs,  or  systems  people  can  help  the  user  mod- 
ify his  program.  Therefore,  the  model  is  written  in  ANSI  standard  Fortran. 


4.  GENERAL  DESCRIPTION  OF  THE  MODEL 
4.1.  Outline  of  the  Model 

The  computer  simulation  model  developed  here  is  a collection  of  computer 
subprograms,  each  of  which,  except  the  one  for  Fourier  transforms,  either 
simulates  a basic  component  of  communication  systems  or  calculates  character- 
istics of  a component.  The  model  does  not  include  a main  program;  a user 
must  write  his  own  main  program,  usually  in  the  Fortran  language,  to  perform 
computer  simulations. 

The  model  includes  19  subprograms;  i.e.,  five  for  modulators,  five  for 
demodulators,  two  for  signal  generators,  two  for  a bandpass  filter,  two  for 
psophometers,  two  for  radio-frequency  (RF)  combining  circuits,  and  one  for 
discrete  Fourier  transforms.  These  19  subprograms  are  listed  in  table  1 with 
respective  brief  descriptions  of  action.  Complete  Fortran  listings  of  the 
subprograms  are  given  in  Appendix  A,  and  their  user  write-ups  are  given  in 
Appendix  B. 
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Table  1.  List  of  Subroutines  Included  in  the  Model 


BPFCH  Simulates  a bandpass  filter  of  Chebyshev  type. 


CBCCT2  Simulates  an  RF  combining  circuit  that  linearly  combines  two 

RF  signals. 

CBCCT3  Simulates  an  RF  combining  circuit  that  linearly  combines  three 

RF  signals. 

DEMAM  Simulates  an  AM  demodulator. 

DEMFM  Simulates  an  FM  demodulator. 

DEMFSK  Simulates  an  FSK  demodulator. 

DEMPHM  Simulates  a phase-modulation  demodulator. 

DEMSSB  Simulates  an  SSB  demodulator. 

DFT  Performs  discrete  Fourier  transform.  1 

FRBPFC  Calculates  frequency  response  of  a bandpass  filter  of 

Chebyshev  type. 

MODAM  Simulates  an  AM  modulator. 

MODFM  Simulates  an  FM  modulator. 

MODFSK  Simulates  an  FSK  modulator.  j 

MOOPHM  Simulates  a phase-modulation  modulator.  j 

MODSSB  Simulates  an  SSB  modulator. 

PSPHCT  Simulates  a psophometer  for  conmercial -telephone  circuit. 

PSPHPT  Simulates  a psophometer  for  program-transmission  circuit. 

SGPLS  Simulates  a pulse-signal  generator. 

SGWGN  Simulates  a white-Gaussian-noise  generator. 


All  subprograms  included  in  the  model  are  written  in  ANSI  standard  For- 
tran (ANSI,  1966).  They  can  be  called  by  a Fortran  program. 

All  signals  and  noise  are  represented  as  their  complex  voltage  spectra 
at  the  input  and  output  of  each  subroutine.  This  means  that,  insofar  as  the 
user  is  concerned,  the  frequency-domain  approach  is  adopted,  although  signals 
are  transformed  to  time  functions  in  some  subroutines.  The  spectrum  of  each 
signal  or  noise  is  a complex  voltage  spectrum. 

The  power  spectral  density  of  the  output  of  the  Gaussian  noise  generator 
subroutine  is  unity.  The  output  powers  of  all  signal  generator  and  modulator 
subroutines  are  normalized  in  such  a way  that  the  output  power  is  equal  to 
unity,  where  the  ITU  Radio  Regulations  (ITU,  1968)  means  for  specifying  the 
power  of  an  emission  are  used. 

i 
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Becaust  o!  the  use  of  a simplified  fast  Fourier  transform  (Cooley  and 
Tukey,  1965;  Cochran  et  al.,  1967),  the  number  of  complex  spectrum  components 
of  the  input  and  output  signals  must  be  a power  of  2.  In  addition,  we  have 
imposed  a restriction  that  the  number  of  complex  spectrum  components  must  not 
exceed  2^*'. 

Some  subroutines  included  in  the  model  call  the  RANF(X)  function,  which 
generates  a pseudo-random  number  uniformly  distributed  in  the  interval  be- 
tween 0 and  1.  (The  parameter  X is  a dummy  parameter.)  When  a user's  com- 
puter system  does  not  include  the  RANF  function  or  its  equivalent  as  a li- 
brary function,  the  user  must  supply  the  function. 

A prospective  user  of  this  computer  simulation  model  must  have  a good 
knowledge  of  communication  systems;  otherwise  he  cannot  perform  meaningful 
simulation.  By  a good  knowledge  of  communication  systems,  we  do  not  mean 
details  of  communication  systems.  What  is  meant  by  a good  knowledge  here  is 
a good  common-sense  knowledge  regarding  the  power  of  emission  (or  transmitter 
power),  noise  power,  noise  density,  signal-to-interference  ratio,  signal-to- 
noise  ratio,  received  signal  quality,  etc.,  as  well  as  good  understanding  of 
meanings  of  system  parameters  for  the  systems  he  plans  to  simulate. 

A prospective  user  of  this  model  must  know  how  to  use  a computer  and 
must  have  a good  knowledge  of  Fortran,  especially  the  ANSI  standard  version 
of  Fortran.  This,  however,  does  not  mean  that  a prospective  user  must  be  an 
expert  in  computer  programming. 

The  size  of  the  computer  required  for  simulation  with  this  model  depends 
on  the  size  of  simulation,  i.e.,  the  number  of  frequency  components  (or 
equivalently,  the  number  of  sampling  points)  for  analog  signals  and  the  length 
of  the  bit  sequence  in  digital  systems.  In  many  applications,  however,  the 
required  computer  size  is  not  a severe  problem.  Since  the  computer  must  have 
a Fortran  compiler  which  requires  a large  core  size,  we  can  generally  say 
that  a computer  can  run  a program  with  this  model  if  it  can  compile  the  model. 


4.2.  Mathematical  Bases 

Mathematical  bases  for  some  subprograms  included  in  the  model  are  des- 
cribed in  this  section.  Modulators  and  signal  and  noise  generators  are  des- 
cribed first,  followed  by  the  descriptions  of  demodulators,  bandpass  filter 
characteristics,  and  psophometers. 

AM  and  SSB  Modulator:  These  subroutines  generate,  from  the  spectrum 

components  of  the  modulating  signal,  the  sideband  components  in  the  frequency 
doma i n . 

FM  Modulator:  This  subroutine  integrates  the  modulating  signal  with  re- 

spect to  time  and  phase  modulates  the  carrier  with  the  integrated  modulating 
signal  in  the  same  manner  as  in  the  phase-modulation  modulator  (described  be- 
low). Integration  with  respect  to  time  is  equivalently  performed  in  the  fre- 
quency domain  by  dividing  each  frequency  component  of  the  signal  by  ju)=j2Trf. 


'FSK  Modulator:  This  subroutine  generates  an  FSK  signal  in  the  time  do- 

main and  Fourier  transforms  the  signal  to  the  frequency  domain. 

Phase-Modulation  Modulator:  This  subroutine  Fourier  transforms  the  mod- 

ulating signal  from  the  frequency  domain  to  the  time  domain,  phase  modulates 
the  carrier  with  the  transformed  signal  in  the  time  domain,  and  Fourier  trans 
forms  the  modulated  signal  to  the  frequency  domain. 

Pulse-Signal  Generator:  This  subroutine  generates  a pulse  signal  in  the 

time  domain  and  Fourier  transforms  the  generated  signal  to  the  frequency  do- 
main. 


White-Gaussian-Noise  Generator:  This  subroutine  generates,  for  each  com 

plex  component  of  the  noise,  a uniformly  distributed  random  number,  converts 
the  number  to  a Rayleigh  distributed  random  number,  and  multiplies  the  number 
by  cosine  and  sine  of  an  arbitrary  (uniformly  distributed)  phase  angle  to  ob- 
tain the  real  and  imaginary  parts  of  the  spectrum  component.  This  procedure 
is  based  on  the  method  suggested  by  Box  and  Muller  (1958)o 

AM  and  Phase-Modulation  Demodulators:  These  subroutines  Fourier  trans- 

form the  input  signal  from  the  frequency  domain  to  the  time  domain,  extract 
the  amplitude  (or  envelope  voltage)  and  the  phase  of  the  signal,  respectively 
in  the  time  domain,  and  Fourier  transform  the  extracted  information  from  the 
time  domain  to  the  frequency  domain. 

PM  Demodulator:  This  subroutine  calculates  the  instantaneous  frequency 

of  the  input  signal  from  the  spectrum  of  the  input  signal  and  Fourier  trans- 
forms the  instantaneous  frequency  from  the  time  domain  to  the  frequency  do- 
main. The  instantaneous  frequency  is  calculated  as  follows: 

We  assume  that  the  input  signal  is  a narrow-band  signal;  i.e.,  the  input 
signal  can  be  expressed  as  a complex  function  of  time  t by 

e(t)  = A(t)  exp[j0(t)J , 

where  A(t)  and  0(t)  are  real  functions  of  t representing  the  instantaneous 
amplitude  and  phase,  respectively.  Taking  the  natural  logarithms  of  both 
sides  of  this  equation,  we  have 

SLn  e(t)  = in  A(t)  + J0(t). 

Since  A(t)  is  a real  function  of  t,  in  A(t)  is  also  a real  function  of  t. 
Therefore,  0(t)  can  be  expressed  by 

0(t)  = lm{in  e(t)}, 

where  Im{...}  stands  for  "the  imaginary  part  of."  Since  the  instantaneous 
frequency  f (t)  is  the  time  derivative  of  the  instantaneous  phase  divided  by 
Zm,  it  is  expressed  by 

""  277  dt  ^ dt 
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When  e(t)  is  a periodic  signal  with  its  fundamental  frequency  fg  and  is  ex- 
pressed by 

e(t)  =LC|^exp[j2iTkfQt] 
k 

with  complex  constants  Cj^,  its  time  derivative  can  be  expressed  by 

= j 2TTfQ  E kC^exp  [ j 2TTkf ^t] . 
k 

From  these  equations,  therefore,  we  have 

f(t)  = Re{fQEkC|^exp[j2Trkfjjt]/Ec,^exp[j27rkfQt]} 
k k 

as  an  expression  for  the  instantaneous  frequency,  where  Re{...}  stands  for 
"the  real  part  of." 

FSK  Demodulator:  This  subroutine  calculates  the  instantaneous  frequency 

of  the  input  signal  from  the  spectrum  of  the  input  signal  in  the  same  manner 
as  in  the  FM  demodulator  and  estimates  the  keying  bit  sequence.  Bit  synchro- 
nization is  established  so  as  to  maximize  the  crosscorrelation  between  the 
sequence  of  instantaneous  frequency  values  and  the  estimated  bit  sequence. 

SSB  Demodulator:  This  subroutine  demodulates  the  input  signal  by  shift- 

ing down  the  spectrum  components  in  the  frequency  domain. 

Calaulation  of  Frequency  Response  of  a Chebyshev-Type  Bandpass  Filter: 

It  is  based  on  the  description  given  by  Storer  (1957).  Necessary  equations 
are  listed  as  follows: 

Frequency  response  of  a Chebyshev-type  bandpass  filter  is  uniquely  spec- 
ified by  the  following  parameters: 

B = 3-dB  bandwidth, 
fc  = center  frequency, 
n = order, 

R = ripple  in  the  passband  in  dB. 

The  complex  frequency  response  F(f)  can  be  expressed  by 

F(f)  = F(fc)  • n 

k=l 

where  F(fc)  = 1 when  n is  odd, 

= 1/Vr  when  n is  even, 

= ripple  in  the  passband  in  voltage  ratio 
= ioR/20. 

Pk  = pole  of  the  filter  (k  = l,2,...,n), 
r = scale  factor. 


The  pole  P|^  on  a complex  plane  is  given  by 

P|^  = -sinh  Uo  cos  4>|^  + j cosh  Uq  sin  4)|^, 


where  uo  = [sinh"Ml/e)]/n, 

(})|^  = (2k-n-l)T7/(2ii). 

The  scale  factor  r is  the  largest  real  root  of 

Tn(r)  = 1/e, 

XU 

where  Tn(x)  is  a Chebyshev  polynomial  of  the  first  kind  of  the  n^"  order 
(Abramowitz  and  Stegjn,  1964) „ This  polynomial  is  defined  by 

Tn(x)  = cos(n  cos"^x). 

It  can  also  be  calculated  from  the  following  recurrence  relation 

Tn(x)  = 2xTn-i(x)  - Tn-2(x) 

with  To(x)  = 1, 

Ti(x)  = X. 

Numerical  solution  is  required  for  determining  the  scale  factor  r for  a given 
value  of  e when  n equals  3 or  greater. 

Psophometers:  These  subroutines  are  based  on  the  weighting-coefficient 

curves,  recommended  by  the  CCITT  (1973)  and  shown  in  figure  2.  Thirty  val- 
ues of  the  weighting  coefficient  at  30  frequencies  are  stored  in  each  subrou- 
tine. When  it  is  called,  each  subroutine  calculates  the  weighting  coeffi- 
cient for  the  frequency  of  each  component  of  the  input  signal  (by  linear 
interpolation  with  respect  to  the  logarithm  of  the  frequency),  multiplies  the 
coefficients  by  the  input  signal  spectrum  on  a compcnent-by-component  basis, 
and  calculates  the  total  output  power  (by  numerical  integration). 


4.3o  Validation 

All  validation  tests  have  been  performed  on  the  CDC-6600  computer  of  the 
U.S.  Department  of  Commerce  Boulder  Laboratories. 

As  validation  tests  of  the  modulators  and  demodulators  of  the  AM,  FM, 
phase-modulation,  and  SSB  systems,  the  back-to-back  tests  for  the  communica- 
tion systems  are  simulated,  in  which  the  modulator  output  signal  is  fed  di- 
rectly to  the  demodulator  input.  An  arbitrary  modulating  signal  is  applied 
to  the  modulator  input,  and  after  modulation  and  demodulation,  the  demodula- 
tor output  signal  (demodulated  signal)  is  compared  with  the  original  modulat- 
ing signal.  The  comparisons  indicate  that  the  corresponding  spectrum  compo- 
nents of  both  signals  agree  in  at  least  5 to  6 decimal  digits. 
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Figure  2.  Frequency  responses  of  the  CCITT  psophometers  for 
commercial -telephone  and  program-transmission  cir- 
cuits. 


( 

i Validation  tests  for  the  modulator  and  demodulator  of  the  FSK  system  are 

I performed  also  as  the  back-to-back  tests.  The  test  result  indicates  that  the 

[ estimated  bit  sequence  from  the  demodulator  always  agrees  with  the  input  bit 

I sequence  given  to  the  modulator  as  long  as  the  number  of  bits  in  the  total 

I time  period  does  not  exceed  one  half  of  the  number  of  frequency  spectrum  com- 

I ponents  of  the  modulated  signal. 

[ 

I Another  validation  test  is  performed  on  the  FM  modulator.  A sinusoidal 

signal  is  given  to  the  FM  modulator  input  as  the  modulating  signal,  and  the 
output  spectrum  is  compared  with  the  FM  spectrum  calculated  theoretically 
(Schwartz,  1970),  Both  spectra  mutually  agree  very  well,  in  at  least  10 
decimal  digits. 

The  Fourier-transform  subprogram  is  also  tested.  The  spectrum  of  an  ar- 
bitrary signal  is  transformed  to  the  time  domain  first,  and  the  time-domain 
signal  is  transformed  back  to  the  frequency  domain.  The  test  result  indi- 
cates that  each  component  of  the  output  spectrum  agrees  with  the  correspond- 
ing component  of  the  input  spectrum  in  10  decimal  digits  or  more. 


<1 
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5o  DW  TO  USE  THE  MODEL 
5.1.  General  Guidelines  for  the  User 

Since  the  model  consists  of  subroutine  subprograms  but  does  not  include 
a main  program,  each  user  must  write  a main  program  that  properly  calls  neces- 
I sary  subroutines.  The  main  program  is  usually  written  in  Fortran.  Since  com- 

I puter  operators  will  run  the  program  on  tne  computer,  writing  ’the  main  program 

is  a major  part  of  computer  simulation  for  the  user  of  the  model.  Although  no 
I • standardized  procedure  exists  for  writing  a main  program,  we  will  describe  a 

I typical  procedure  that  can  easily  be  followed  by  prospective  users. 

i 

j • In  many  cases,  the  user  must,  as  the  first  step,  select  necessary  sub- 

i routines  for  the  desired-signal  source  and  the  interfering-signal  and/or  noise 

sources.  Selection  of  these  subroutines  depends  on  the  problem  the  user 
wishes  to  solve  or  the  system  he  wishes  to  simulate.  Any  subroutine  that  sim- 
ulates a nwdulator  or  a pulse  signal  generator  can  be  a desired-signal  source. 
Any  subroutine  that  simulates  a modulator  or  a signal  or  noise  generator  can 
be  an  interfering-signal  or  a noise  source.  A subroutine  can  be  used  as  a 
desired-signal  source  and  an  interfering-signal  source  at  the  same  time.  More 
than  one  subroutine  can  be  used  simultaneously  as  interfering-signal  sources. 
The  output  signals  from  these  subroutines  are  given  as  complex  voltage  spec- 
tra. As  described  earlier,  the  power  spectrum  density  of  the  noise  generator 
output  is  unity,  and  the  output  powers  of  all  other  subroutines,  expressed  in 
terms  of  what  is  used  in  specifying  the  power  of  an  emission  in  the  ITU  Radio 
Regulations  (1968),  are  equal  to  unity,  this  will  make  the  adjustment  of  the 
signal -to-interference  ratio  very  easy  and  straightforward.  Although  these 
subroutines  have  a number  of  formal  parameters,  it  is  necessary  only  to  iden- 
tify each  parameter;  no  dimensions  of  arrays  or  numerical  values  of  variables 
have  to  be  determined  at  this  stage. 

The  next  step  is  to  simulate  the  receiver  input  or  to  combine  the  inter- 
fering signal  and/or  noise  with  the  desired  signal.  The  model  includes  two 
subroutines  for  RF  signal  combiners;  one  linearly  combines  two  signals  with  a 
specified  ratio,  and  the  other  linearly  combines  three  signals  with  two  speci- 
fied ratios.  Again,  no  numerical  values  have  to  be  assigned  to  these  ratios 
at  this  stage. 

The  next  step  is  to  simulate  radio-frequency  (RF)  and  intermediate- 
frequency  (IF)  amplifiers  with  a bandpass  filter.  The  model  includes  a sub- 
routine that  simulates  a Chebyshev-type  bandpass  filter.  When  this  subrou- 
tine is  called  with  new  values  of  parameters,  it  calls  another  subroutine 
that  is  also  included  in  the  model  and  calculates  the  frequency  response  of 
the  bandpass  filter. 

The  next  step  is  to  demodulate  the  bandpass  filter  output  signal  with  a 
demodulator  subroutine  that  corresponds  to  the  desired  signal.  The  output 
from  the  FSK  demodulator  is  a sequence  of  bits,  while  the  outputs  from  all 
other  demodulator:  are  given  as  complex  voltage  spectra  in  the  baseband. 

From  the  voltage  spectrum  in  the  baseband,  one  can  derive  a measure  of 
quality  of  the  output  signal.  The  model  includes  two  subroutines  that  simu- 
late the  psophometers  for  the  commercial -telephone  and  program- transmission 
circuits,  recommended  by  the  CCITT  (1973). 
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The  steps  described  in  the  above  five  paragraphs  constitute  a main  part  I 

of  the  program.  This  part  can  be  repeated  with  various  parameter  values  in  a | 

DO-loop  or  DO-loops.  In  some  applications,  this  part  must  be  executed  twice,  i 

with  and  without  the  modulating  signal. 

I 

Each  subroutine  for  a modulator,  except  the  one  for  the  FSK  modulator,  s 

accepts  the  input  modulating  signal  in  the  form  of  a voltage  spectrum.  Also, 
each  modulator  or  signal  generator  gives  its  output  in  the  form  of  a voltage 
spectrum.  The  user  must  select  suitable  values  for  the  number  of  frequency 
components,  NFC,  and  the  unit  frequency  spacing,  UFS,  in  the  spectrum.  Selec- 
tion of  proper  values  for  NFC  and  UFS  requires  some  engineering  judgments. 

The  following  descriptions  may  help  the  user  make  such  selections.  Each 
signal  appearing  in  the  model  is  a periodic  function  of  time;  the  period  of 
the  function  is  represented  by  the  reciprocal  of  UFS.  The  frequency  of  the 
highest  frequency  component  that  can  appear  in  the  model  is  (NFC-1)*UFS.  The 
unit  time  interval  of  the  sampling  points  is  represented  by  1/(NFC*UFS).  When 
harmonic  distortion  in  the  demodulator  output  signal  is  discussed,  selection 
of  NFC*UFS  i 2*FRMX  is  recommended,  where  FRMX  is  the  frequency  at  which  the 
harmonic  distortion  component  essentially  vanishes. 

In  order  to  save  computation  time  required  for  the  Fourier  transform, 
the  model  imposes  a restriction  on  the  value  of  NFC  that  NFC  be  equal  to  a 
power  of  two. 
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5.2.  An  Example  — AM  Signal  With  AM  Interference 

^ Let  us  illustrate  the  use  of  the  model  with  a simple  example.  In  this 

j example,  we  will  determine  the  output  signal-to-interference  ratio  (SIR)  when 

1 an  amplitude-modulated  (AM)  desired  signal  is  interfered  with  by  an  AM  signal 

in  the  absence  of  noise.  The  output  SIR  depends  on  the  modulation  indexes  of 
both  the  desired  and  interfering  signals,  the  frequency  difference  between 
the  carriers  of  the  two  signals,  and  the  radio-frequency  (RF)  SIR.  For  sim- 
plicity we  do  not  include  the  effect  of  a bandpass  filter  in  this  example. 

We  measure  the  powers  at  the  commercial -telephone  psophometer  output  with 
and  without  modulation  of  the  desired  signal  and  calculate  the  output  SIR  as 
the  ratio  of  the  difference  between  the  two  output  powers  to  the  latter  (i.e., 
without  modulation). 

! The  first  step  it  to  review  briefly  what  subroutines  are  needed.  Since 

1 both  the  desired  signal  and  the  interfering  signal  are  AM  signals,  we  need 

I only  the  MODAM  subroutine  as  the  modulator.  We  can  use  the  CBCCT2  subrou- 

tine to  superpose  the  interfering  signal  upon  the  desired  signal  because  we 
have  only  two  signals  to  combine  in  this  example.  We  use  the  DEMAM  subrou- 
tine to  demodulate  the  combined  signal.  To  calculate  the  psophometrical ly- 
weighted  signal  power  from  the  demodulator  output  signal,  we  use  the  PSPHCT 
subroutine  that  simulates  the  psophometer  for  a commercial -telephone  circuit. 
Thus,  the  model  includes  all  necessary  suoroutines  for  this  example. 

We  will  use  numerals  1 and  2 for  the  desired  signal  and  the  interfering 
signal,  respectively,  as  the  endings  of  the  variable  and  array  names.  For 
example,  FCl  denotes  the  center  frequency  of  the  desired  signal,  and  SF2  is 
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the  name  of  the  array  for  the  signal  that  modulates  the  interfering  signal. 
Then  the  following  sequence  of  Fortran  call  statements  is  almost  self-evident: 

C GENERATES  THE  MODULATED  DESIRED  SIGNAL. 

CALL  MODAM  (NFC,UFS,SFl,FCl,MIi,VFl) 

C GENERATES  THE  INTERFERING  SIGNAL  (MODULATED). 

CALL  MODAM  (NFC,UFS,SF2,FC2,MI2,VF2) 

C COMBINES  THE  MODULATED  DESIRED'  SIGNAL  WITH  THE 
C INTERFERING  SIGNAL. 

CALL  CBCCT2  (NFC,VF1,VF2,SIR,VF4) 

C DEMODULATES  THE  COMPOSITE  SIGNAL. 

CALL  DEMAM  (NFC,VF4,SF4) 

C APPLIES  THE  DEMODULATED  SIGNAL  TO  THE  PSOPHOMETER 
C INPUT. 

CALL  PSPHCT  (NFC,UFS,SF4,PS4) 

(For  the  interpretation  of  parameters  for  these  subroutines,  see  Appendix  B.) 
We  do  not  have  to  assign  any  numbers  to  the  input  parameters  or  determine  the 
dimensions  of  arrays  at  this  stage.  In  this  example,  we  must  also  calculate 
the  psophometer  output  power  for  an  unmodulated  desired  signal.  Therefore, 
v;e  add  to  the  above  sequence  of  call  -statements  the  following  statements: 

C GENERATES  THE  UNMODULATED  DESIRED  SIGNAL. 

CALL  MODAM  (NFC,UFS,SF1,FC1,0.0,VF0) 

C COMBINES  THE  UNMODULATED  DESIRED  SIGNAL  WITH  THE 
C INTERFERING  SIGNAL. 

CALL  CBCCT2  (NFC,VF0,VF2,SIRrVF3 ) 

C DEMODULATES  THE  COMPOSITE  SIGNAL. 

CALL  DEMAM  (NFC ,VF3, SF3 ) 

C APPLIES  THE  DEMODULATED  SIGNAL  TO  THE  PSOPHOMETER 
C INPUT. 

CALL  PSPHCT  (NFC,UFS,SF3,PS3) 

(Note  that  generation  of  the  interfering  signal  is  not  duplicated  in  this  sec- 
ond sequence.)  The  order  of  the  above  nine  call  statements  can  be  rearranged 
as  long  as  their  logical  order  is  not  violated;  the  only  restriction  is  that, 
when  the  output  of  a subroutine  A is  the  input  of  another  subroutine  B,  the 
call  to  A must  precede  the  call  to  B. 

From  the  psophometer  output  powers  with  and  without  modulation,  PS4  and 
PS3,  we  can  calculate  the  output  SIR,  SIRO,  as  (PS4-PS3)/PS3  or  as  PS4/PS3-1. 
To  avoid  a possible  overflow,  we  set  SIRO  = 10^°  (or  200  dB)  when  PS3  is  zero. 
To  avoid  calling  the  ALOG10  function  (cotimon  logarithmic  function)  with  a non- 
positive argument  in  converting  SIRO  -from  a power  ratio  to  a dB  value  we  can 
modify  the  SIRO  value  when  it  is  less  than  a certain  value;  we  set  this  value 
to  be  10"^“  (or  -100  dB).  These  two  things  are  minor  technical  points  but 
are  sometimes  serious  enough  to  frustrate  the  programmer.  Finally,  we  con- 
vert SIRO  from  a power  ratio  to  a dB  value.  The  procedure  described  in  this 
paragraph  can  be  summarized  in  the  following  sequence  of  Fortran  statements: 

IF(PS3.LE.0.0)  GO  TO  11 
SIRO  = PS4/PS3-1.0 
IF(SIr6.LT.1.0E-10)  SIRO  = 1.0E-10 
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GO  TO  12 

11  SIRO  = 1.0E+20_ 

12  SIRO  = 10.0*ALOG10(SIRO) 

This  sequence  of  the  Fortran  statements  given  in  the  preceding  two  para- 
graphs is  the  core  of  a program  that  simulates  reception  of  an  AM  signal  in- 
terfered with  by  another  AM  signal.  Except  for  some  changes  in  the  order  of 
the  statements,  this  core  is  common  in  many  such  programs.  This  sequence  can 
be  placed  in  a DO-loop  or  DO-loops  when  simulation  is  done  for  more  than  one 
set  of  parameter  values. 

What  a user  of  the  model  does  next  with  this  core  portion  of  the  program 
depends  on  his  purpose.  A user  may  be  interested  in  the  effect  of  frequency 
difference  between  the  desired  and  interfering  signals,  while  another  user 
may  be  interested  in  only  one  value  of  frequency  difference.  A user  may  wish 
to  write  the  main  program  in  a closed  form  with  all  input  data  included  in 
the  program,  while  another  user  may  write  the  program  in  an  open  form  so  that 
the  program  will  read  in  necessary  data  during  execution. 

Table  2 lists  the  AMAM00  program  that  is  written  as  an  example  for  illus- 
tration purposes.  This  program  simulates  reception  of  an  AH  desired  signal  in- 
terfered with  by  another  AM  signal  for  two  values  of  modulation  index  for  the 
desired  signal  (fill  = 0.5,  1.0),  for  two  values  of  modulation  index  for  the  in- 
terfering signal  (MI2  = 0.5,  1.0),  for  31  values  of  frequency  difference  be- 
tween the  desired  and  interfering  signals,  and  for  seven  values  of  the  input 
RF  SIR.  These  combinations  of  parameters  are  achieved  by  programming  with 
nested  4-deep  oO-loops.  This  program  is  a closed-form  program  and  includes 
all  necessary  data  in  it.  Values  of  the  frequency  difference  (OF)  and  of  the 
RF  SIR  are  calculated  in  the  beginning  part  of  the  program,  while  all  other 
data  are  given  to  the  program  with  the  data  initialization  statements.  This 
program  assumes  that  the  desired  and  interfering  signals  are  modulated  with 
1000-Hz  and  1200-Hz  tones,  respectively.  We  have  selected 

NFC  = 128, 

UFS  = 200  (Hz), 

FCl  = 10000  (Hz), 

where  NFC  is  the  number  of  frequency  components  in  the  spectrum  of  a signal, 

UFS  is  the  unit  frequency  spacing  in  the  spectrum,  and  FCl  is  the  center  fre- 
quency of  the  desir^  signal.  Use  of  this  UFS  value  limits  all  frequency  val- 
ues used  in  this  program  to  integral  multiples  of  200  Hz.  Use  of  this  NFC 
value  assures  that  all  nonzero  components  of  the  modulated  signals  are  in- 
cluded in  the  VF  arrays,  because  the  VF  arrays  cover  the  frequency  range  from 
0 to  200  X (128-1)  = 25400  Hz. 

A part  of  the  result  from  the  AMAM00  program  is  shown  in  table  3.  The 
AMAM00  program  produces  a total  of  four  pages  of  computer  printouts  corres- 
ponding to  the  four  combinations  of  two  modulation  index  values  for  the  de- 
sired signal  and  two  modulation  index  values  for  the  interfering  signal.  The 
computation  time  required  to  run  the  AMAM00  program  on  the  CDC-6600  computer 
was  approximately  60  seconds  (at  a cost  of  approximately  $8).  Tafe  3 cor- 
responds to  the  case  where  both  indexes  are  equal  to  unity.  In  this  table, 

Mil  and  MI2  are  the  modulation  indexes  for  the  desired  signal  and  interfering 
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Table  2.  Fortran  Listing  of  the  AMAM00  Program 


PROGRAM  AMAMCC (OUTPUT, TAPE6=OUTPUT) 

C Ah'  SYSTEM  WITH  AM  IMERFERENCE  AND  NO  NOISE 
C THE  MODULATING  SIGNAL  FOR  THE  CESIREO  SIGNAL  IS  A 1003-HZ 
C SINLSOIO,  AND  THAT  FOR  THE  INTERFERING  SIGNAL  IS  A 1200-hZ 
C SINUSOID. 

c declaration  statements 

REAL  Mil, MI2, Mill, MI2I 

DIMENSION  MIK2)  ,MI2(  2)  ,DF(31)  ,SIR(7)  ,SIfiO(D 
COMPLEX  SF1(128),SF2(128),SF3(128),SF4(128), 

1 VFO (128) ,VF1(128)  ,VF2(128I  ,VF3 (128», VF4(126I 

* DATA  NMl/2/,  NM2/2/,  NDF/31/,  NSIR/7/ 

DATA  NFC/128/,  UFS/200,0/,  FCl /1 0 0 0 0 . 0 / , 

1 MI1/0.F,1.0/,  MI2/0.5,1.0/, 

2 SF1/5*((0.0,0.0)),(1.0,0.0),122*((0. 0,0.0))/, 

3 SF2/G«((0.0,0.0)),(1.0,0.0),121«((0.0,0.0))/ 

DATA  NAME/6HAMAM00/,  LUN/6/ 

C DATA  PREPARATION 

10  DO  11  IDF=1,NCF 
OF(IOF)=200* (IOF-1) 

IF( IDF.GT.2F)  DF( IOF)=1000»(IDF-21) 

11  CONTINUE 
DO  12  1SIR=1,NSIR 

SIR (ISIR)  = 5*  IISIR-1) 

12  CONTINUE 
C SIMULATION 

20  DO  89  IM1=1,NM1 

C - GENERATES  THE  UNMOCULATED  DESIRED  SIGNAL, 

CALL  MODAM(NFC,UFS,SF1,FC1,0.0,  VFO) 

C - GENERATES  THE  MODULATED  DESIRED  SIGNAL. 

MI1I=MI1 (IMl) 

CALL  MOOAM(NFC»UFS,SFl,FCl,MIlI, VFl) 

DO  83  IM2=1,NM2 

C - GENERATES  THE  INTERFERING  SIGNAL  (MODULATED). 

MI2I=MI2 (IM2) 

DO  79  IOF=l,NOF 
FC2=FC1+CF(IDF) 

CALL  M0CAM(NFC,UFS,SF2,FC2,MI2I  ,VF2) 

DO  59  1SIR=1,NSIR 

C - COMBINES  the  UNMODULATED  DESIRED  SIGNAL  WITH  THE  INTERFERING 

c signal. 

CALL  CECCT2(NFC ,VF0,^F2,SIR( ISIR) , VF3) 

C - DEMODULATES  THE  COMPOSITE  SIGNAL. 

CALL  OEMAM (NFC, VF3  ,SF3) 

C - APPLIES  THE  DEMODULATED  SIGNAL  TO  THE  PSOPHOMETER  INPUT. 

CALL  FSPHCT(NFC ,UFS,SF3,PS3) 

C - CCM6INES  THE  MODULATED  DESIRED  SIGNAL  WITH  Tt-E  INTERFERING 
C SIGNAL. 

call  CBCCT2(NFC, VF1,VF2,SIR( ISIR), VF4) 

C - DEMODULATES  THE  COMPOSITE  SIGNAL. 

CALL  DEMAM(NFC, VF4,SF4) 

C - APPLIES  THE  OEMC'DULATEO  SIGNAL  TO  THE  PSOPHOMETER  INPUT. 

CALL  FSPHCT (NFC ,UFS,SF4,PS4) 

C CALCULATION 

40  IF(PS3.LE. 0. 0)  GC  TO  42 

SIROI=FS4/PS3-1.0 

IF(SIFCI.LE.l.OE-lO)  SIROI= 1 . OE-10 
GO  TO  43 
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Table  3.  Sample  Printout  Produced  by  the  AMAM00  Program, 
Listed  in  Table  2. 

4M4MOO  AM  SYSTEM  WITH  AM  INTEPFERENCE 

MI1  = 1.C0(J  MI2  = 1.00C 


OUTPUT  SIGKAL-TO-INTERFEPENCE  ratio  (03) 


OF 

INPUT  SIR  (03)  = 

0*0  5*0  llj*J 

15.0 

?0.0 

25.0 

30.0 

0. 

2.9 

8.1 

13.5 

19.0 

24.4 

29.7 

34.9 

20  0. 

1.6 

7.6 

13.0 

16.4 

23.6 

28.7 

33.  7 

400. 

-.2 

6.0 

11.7 

17.3 

22.6 

27.7 

32.8 

600. 

-1.2 

5.8 

11.7 

17.2 

22.4 

27.6 

32.6 

eoo. 

-.5 

4.3 

9.6 

14.9 

20.1 

25.2 

30.2 

100  0. 

-100.0 

-2.0 

6.9 

13.3 

18.9 

24,2 

29.4 

1200. 

-100.0 

1.9 

9.0 

14.9 

20.2 

26.4 

30.5 

1400. 

-2.6 

4.6 

10.5 

15.9 

21.2 

26 . 3 

31.4 

1600. 

-1.2 

5.2 

10.9 

16.4 

21.7 

26 .8 

31.8 

1000. 

-.9 

5.1 

11.0 

16.5 

21.8 

26.9 

31.9 

200  0. 

-2.6 

4.3 

10.5 

16.3 

21.6 

26.8 

31.9 

2200. 

1.0 

6,1 

11,4 

16,7 

21,9 

26.9 

32.  0 

2400. 

-1.0 

5.2 

11.3 

16.9 

22.3 

27.5 

32.6 

2600. 

-.1 

6.6 

12.4 

17.9 

23.2 

28.3 

33.4 

2800. 

.3 

7.3 

13.1 

18.6 

23.9 

29.0 

34.1 

3 0 0 0. 

.2 

7.6 

13.6 

19.2 

24.5 

29,6 

34.7 

3200. 

.1 

8.2 

14.5 

20.1 

2E  .5 

30.6 

35.7 

3400. 

.4 

9.2 

15.6 

21.2 

2E  .4 

31.5 

36.6 

3600. 

.2 

10.0 

16.8 

22.4 

27.6 

32.7 

37.8 

3800. 

1.5 

10.  3 

17.3 

2 3.6 

29.0 

34.2 

39.2 

4000. 

1.6 

11.4 

18.8 

24.8 

30.2 

35.4 

40.5 

4200. 

1.7 

11.9 

19.8 

26.0 

31.5 

36.7 

41.8 

4400. 

1.4 

12.5 

21.0 

27.5 

33,1 

36.3 

43.4 

4600. 

1.8 

12.9 

21.9 

28.6 

34.3 

39.6 

44.7 

4800. 

1.8 

13.1 

22.8 

30.0 

36.0 

41.3 

46.  5 

5000. 

1.8 

13,6 

23.8 

31.8 

38,2 

43.8 

49.0 

6000. 

1.7 

13.9 

25.7 

36 . 3 

46.4 

56.0 

64.7 

7000. 

2.0 

14.5 

26.1 

36.5 

46.7 

56.8 

66.9 

8000. 

1.8 

13.8 

26.5 

36.3 

46.7 

56.8 

66.9 

9000. 

1.9 

13.8 

25.5 

36.3 

46.6 

56.8 

66.9 

000  0. 

1.9 

13.6 

25.5 

36.3 

46.6 

56.8 

66.8 

) 
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C i-.Mi.  .LATES  THE  FILTER  OUTPUT  SIGNAL. 

CALL  DEMFM  (NFC,UFS,VF2,FC,SF2) 

(For  the  interpretation  of  parameters  for  these  subroutines,  see  Appendix  B.) 
Again,  no  numbers  have  to  be  assigned  to  any  parameters  at  this  stage. 

The  above  sequence  of  call  statements  leaves  the  spectrum  of  the  demodu- 
lated signal  in  the  SF2  array  as  the  result  of  simulation.  Although  it  may 
look  a little  complicated,  calculation  of  the  final  result  from  the  result  of 
the  simulation  left  in  SF2  is  by  no  means  difficult  once  one  understands  the 
physical  phenomena. 

If  the  modulated  signal  passes  through  no  bandpass  filter,  the  SF2  array 
will  coincide  with  the  SFl  array  that  represents  the  modulating  signal  spec- 
trum. If  the  modulated  signal  passes  through  a bandpass  filter,  on  the  other 
hand,  SF2  generally  differs  from  SFl.  Any  difference  between  SF2  and  SFl  is 
the  result  of  the  effect  of  the  bandpass  filter.  The  way  of  representing  the 
difference  between  the  two  spectra  may  not  be  unique  but,  when  modulation 
is  with  a single  tone  as  assumed  in  this  example,  the  following  representa- 
tion seems  reasonable. 

Since  we  have  assumed  modulation  with  a single  tone,  the  SFl  array  has, 
as  the  only  nonzero  element,  its  second  element  that  corresponds  to  the  fun- 
damental frequency  component.  The  second  element  of  the  demodulated  signal 
spectrum  SF2  is  generally  different  from  that  of  SFl  both  in  amplitude  and 
phase,  i.e.,  it  is  usually  suppressed  in  amplitude  and  delayed  in  phase.  We 
define  t signal -suppression  ratio  and  phase  delay  to  rep>"esent  these  dif- 
ferences. In  addition  to  the  differences  in  the  amplitude  and  phase  of  the 
fundamental -frequency  component,  we  can  define  the  total  distortion,  as  the 
ratio  of  the  power  of  the  fundamental  component  to  the  total  power  of  all  har- 
monic components,  and  also  the  harmonic  distortions,  each  as  the  ratio  of  the 
power  of  the  fundamental -frequency  component  to  the  power  of  each  harmonic 
component. 

Harmonic  distortion  powers  and/or  the  total  distortion  power  can  be  zero. 
To  avoid  calling  the  ALO610  function  (common  logarithmic  function)  with  a 
zero  argument  in  converting  the  signal -to-distortion  ratio  from  a power  ratio 
to  a dB  value,  we  set  the  ratio  to  be  1000  dB  when  the  distortion  power  is 
zero. 


The  FMDF00  program  that  is  programmed  for  illustration  purposes  is  listed 
in  table  4„  It  is  again  a closed-form  program;  all  data  are  either  given  with 
the  data  initialization  statements  or  calculated  during  execution.  We  have 
selected 

NFC  = 128, 

UFS  = loO, 

FC  = 64.0, 

FM  = 1.0, 

where  NFC  is  the  number  of  frequency  components  in  the  spectrum  of  a signal, 
UFS  is  the  unit  frequency  spacing  in  the  spectrum,  FC  is  the  center  frequency 
of  the  modulated  signal,  and  FM  is’ the  frequency  of  the  modulating  tone.  For 
simplicity  we  have  normalized  all  frequencies  and  bandwidths  in  the  program 
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Table  4.  Fortran  Listing  of  the  FMDF0(J  Program 


Pf?0G9H1  FHOFtC  I OUTPUT  , T AP  E€  = OUTPUT  ) 

C Ft>  DISTORTION  CAU'FO  BY  A BANDPASS  FILTER  OF  CHEBYSHEw  TYPE 
C DECLARATION  STATEMENTS 

COMPLEX  SFlU2<<>,VFH12‘ntVF?tl2e)  ,SF2  «129) 

DIMENSION  PFtl28»,DSIlO) 

DIMENSION  HI<4)  ,NOB(2  ).»RPE  (3)  ,CF0H  ( 7)  ,CFDF  (4) 

REAL  MI 

DATA  SFl/ (0.0. 0.0 1,(1. Of  0.0) ,126*  <(0.0, 0.0))/ 

DATA  NFC/12F/,  UFS/l.O/t  FC/64.0/,  FM/1.0/, 

1 MI/1.0.2.0,5.0,10.0/. 

2 NOB  /2,‘-/, 

3 RPB  /O. 0,1. 0,2.0/, 

CF3W/0.B,1.0,1.2,l.<<,l.b,l.l,2.0/, 

5 CFOF/O.C ,0.1,0.2, 0.3/ 

DATA  NAME/6HFMOFOO/,  LUN/6/ 

C.  SIMULATION 

10  00  89  IMI=1,»* 

C - GENERATES  AM  FM  SIGNAL. 

F0V  = FM*MI(IMH 
PO=FOV*FOV/2.0 

CALL  MOtFM(KFC,UFS,SFl,FC,FDV,VFl) 

C - DETERMINES  PARAMETERS  FOR • A BANDPASS  FILTER. 

BHCS=2. C» (FM+FCV) 

DO  88  INOB=l,2 
NOBO=NOB(INOE) 

DO  87  IRPB=1,3 
RPB0=RPB(IRP3) 

IF(NOBO.EQ.l. AND. RPBO.NE. 0.01  GO  TO  87 

DO  E9  IBHB=1,7 

8HO=BWCS*CFBH(I0W3) 

DO  59  IFCB=1,4 

OFC=BWB*CFDF( IFCB) 

FCB=FC+OFC 

C - APPLIES  THE  FM  SIGNAL  TO  THE  BANDPASS  FILTER  INPUT. 

CALL  BPFCHtNFC.UPS.NOBO.RPBO ,FCB,BH3, VFl, VF2» 

C - DEMODULATES  THE  FILTER  OUTPUT  SIGNAL. 

CALL  OEMFM(NFC,UFS, VF2,FC,SF2) 

C CALCULATION 

C - CALCULATES  THE  FUNDAMENTAL  FREQUENCY  COMPONENT  OF  THE  OUTPUT 
C SIGNAL  POKER. 

30  P1=(PEAL(SF2(2) > ♦*2*AIMAG ( SF2 < 2) )**2>/2.0 

C - CALCULATES  THE  HARMONIC  COMPONENTS  OF  THE  OUTPUT  SIGNAL 
C POKER  AND  THE  TOTAL  DISTORTION  POKER. 

PT=C.O 

DO  31  1=3, NFC 

PF(I»  = (REAL(SF2{in**2  + AIMAG(SF2<in**2)/2.0 
PT  = PT-fPF(I) 

31  CONTltUE 

C - CALCULATES  THE  SIGNAL  SUPPRESSION  RATIO. 

SSR=-10. 0*ALOG10 (Pl/PO) 

C - CALCULATES  THE  PHASE  DELAY. 

PHS=ATAN2( AIMAG (SF2 (2»» .REAL  (SF2 (2)  ) I 
C - CALCULATES  THE  SIGNAL-TO-OISTCRTION  RATIC. 

IF( PT.GT.O.O)  DT=-10.0*ALOG10(PT/P1) 

IF(FT.LE.O.O)  OT=1000,0 

DO  32  1=3,6 

IF(PF(I) .GT.0.0)  OS(I-l)= 
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Table  4.  Continued 


-10, 0*ALC&10 (eF(  I)  /Fl) 
IF(PF(IJ .LE.O.O)  DS(I-1 ) =10C0.0 
32  CONTINUE 

C PRIMING  OF  THE  RESULTS 

<•0  IFMBHe.EQ.l.  ANO.IFCB.EQ.l) 

1 WRITE  (LUN,2040)  NAME  , M I <1MH  , NC:?0 , K.FBO 

IF(IFCe.GT.l)  GC  TO  42 

41  WRITE  (LUN,2041)  BMB, OFC , SSR , FHS , DT , 

1 <0S(IH) ,IH=2,5) 


GO  1C  59 

42  WRITE  (LUN,2042)  CFG.  SSR,  P HS  , OT  , ( OS  ( IH  ) , IH  = 2 ,5‘) 

59  CONTINUE 

69  CONTINUE 

87  continue 

88  CONTINUE 

89  continue 
CALL  EXIT 

C FORMAT  STATEMENTS 

2040  FORMATdHl.AE  ,8X,36HFM  DISTORTION  DUE  TO  BANDPASS  FUTEF// 

1 15X,4HMI  =,F5.1 ,6X,5HNO0  =,I2,6X,9HRPB(0BI  =,F4.1//// 

2 21X,3HSSR,2X,5HPHASE ,6X,2HOT,5X,2HD2,5X,2HD3, 

3 SX,2H04,5X,2H05/ 

4 3X,3HBH0,5X,3HOFC/ 

5 20X,4H (OB> ,2X,5H(RAO) ,2X,5(3X,4H(0B)>) 

2041  FORMAT<lX/lX,Fe.3,F8.3,2X ,2F7.2,2X,5F7,1) 

2042  FORMAT (IX, 6X,  , F 8. 3, 2X, 2F 7 . 2 , 2X , 5 F7 . 1 ) 

ENO 


to  FM  by  taking  FM  = l„n;  this  convenient  normalization  is  possible  when  a 
program  does  not  simulate  the  psophometer  that  is  dependent  on  the  absolute 
frequency  values.  This  program  performs  simulation  for  a total  of  24  combi- 
nations of  four  values  of  modulation  index  (MI  = 1.0,  2.0,  5.0,  10.0),  two 
values  of  order  of  the  bandpass  filter  (NOB  = 2,  4),  and  three  values  of  rip- 
ple in  the  filter  passband  (RPF  = 0.0,  1.0,  2.0  dB).  For  each  combination, 
FMDF00  calculates  seven  values  of  filter  bandwidth,  BWB,  i.e.,  90%,  100%, 

120%,  140%,  160%,  180%,  and  200%  of  the  Carson  bandwidth,  2*(MI  + 1.0)*FM. 

For  each  BWB  value,  FMDF00  further  calculates  four  values  of  detuning  fre- 
quency, OFC,  as  0%,  10%,  20%,  and  30%  of  the  BWB  value.  (Detuning  frequency 
is  the  difference  in  the  center  frequencies  of  the  bandpass  filter  and  the 
signal,)  For  each  combination  of  MI,  nOb,  and  RPB  and  for  each  value  of  BWB 
and  OFC,  the  FMDF00  program  performs  simulation  and  calculates  the  signal - 
suppression  ratio  (SSR),  phase  (PHASE),  total  distortion  (DT),  and  the  second 
through  the  fifth  harmonic  distortions  (D2,  D3,  D4,  D5). 

A part  of  the  result  from  the  FMDF00  program  is  shown  in  table  5,  This 
table  represents  the  result  for  the  modulation  index  of  ten  (MI  = 10,0),  the 
fourth-order  bandpass  filter  (nOb  = 4),  and  no  ripple  in  the  filter  response 
(RPB  = 0,0  dB),  The  FMDF00  program  yields  23  more  pages  of  similar  printouts 
for  different  combinations.  The  computation  time  required  for  the  FMDF00  pro- 
gram on  the  CDC-6600  computer  was  approximately  48  seconds. 

Table  5 shows  the  signal-suppression  ratio  (SSR),  the  phase  (PHASE),  the 
total  distortion  (DT),  and  the  four  harmonic  distortions  (D2,  03,  D4,  and  D5), 
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Table  5»  Sample  Printout  Produced  by  the  FMDF00  Program,  Listed  in  Table  4. 


FMOF  00 


FM  DISTORTION  DUE  TO  BANDPASS  FILTER 


MI 

= 10.0 

NOB 

= 4 

RPB(CO) 

= 0. 

0 

SSR 

PHASE 

DT 

02 

C3 

’ 04 

05 

BWB 

OFC 

(DB) 

(RAO) 

(OB) 

(OB) 

(OB) 

(03) 

(OB) 

17.600 

0.000 

.08 

-.35 

26 . 1 

2b8. 8 

27.6 

270.6 

35.1 

1.760 

.06 

-.35 

25.6 

36.6 

31.1 

30.7 

33.9 

3.520 

.02 

-.3<« 

21.1 

36.  3 

30.9 

24.6 

31.0 

5.280 

-.02 

-.32 

15.3 

32.6 

24.0 

24.3 

27.3 

72.000 

3.000 

.03 

-.27 

29.  ( 

272.2 

30.1 

271.5 

42.0 

2.200 

.04 

-.27 

2 8.  » 

32.8 

32.3 

41.4 

45.1 

4.400 

.03 

-.27 

27.4 

31.1 

40.6 

32.4 

35.9 

6.600 

.01 

-.27 

23 . 7 

37.5 

28.5 

28.4 

37.5 

26.<»Q0 

0.300 

.01 

-.21 

34.2 

271.7 

34.4 

270.5 

47.9 

2.640 

.02 

-.22 

31.7 

34.1 

35.2 

47.3 

51.6 

5.230 

.03 

-.23  ■ 

29.4 

30.6 

42.7 

40.0 

40.9 

7.920 

.02 

-.24 

28. C 

34.3 

33.6 

32.2 

42.8 

30.800 

0.000 

.00 

-.18 

38.  R 

272.6 

38.9 

285.8 

54.3 

3.080 

.01 

-.18 

34. 2 

36.6 

38.4 

50.6 

56.4 

6.160 

.02 

-.20 

30.9 

31.6 

41.4 

46.8 

46.2 

9.240 

.02 

-.21 

30.4 

33.2 

39.4 

35.9 

44.9 

35.200 

0.000 

.00 

-.15 

43.0 

270.7 

43.1 

270.8 

60.6 

3.520 

.01 

-.16 

37.1 

39.2 

41.6 

54.1 

61.0 

7.040 

.02 

-.17 

32.5 

33.2 

42.0 

52.3 

51.4 

10.560 

.02 

-.18 

32.0 

33.2 

45.9 

39.5 

47.2 

39.600 

0.000 

.00 

-.14 

46.7 

273.1 

46.7 

273.2 

66.5 

3.960 

. 00 

-.14 

39.  8 

41.7 

44.8 

57.7 

65.5 

7.920 

.01 

-.15 

34.2 

34.9 

43.4 

56.2 

56.3 

11.380 

.02 

-.16 

33.1 

33.7 

53.9 

43.0 

49.9 

44.000 

0.000 

.00 

-.12 

49.  9 

271.2 

49.9 

277.1 

71.8 

4.400 

.00 

-.12 

42.3 

43.9 

47.7 

61.3 

69.8 

8 .800 

.01 

-.13 

36.  0 

36.6 

45.2 

59.2 

61.1 

13.200 

.02 

-.15 

34.2 

34.5 

64.5 

46.3 

52.8 

i 


1 

I 
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i for  four  values  of  detuning  frequency  (DFC)  and  seven  values  of  the  filter 

[ bandwidth  (BWB).  Since  the  frequency  of  the  modulating  tone,  FM,  is  taken  to 

f be  unity  in  this  program,  we  can  consider  that  all  frequencies  and  bandwidths 

[ are  normalized  to  the  modulating-tone  frequency.  Although  the  center  fre- 

[ quency  of  the  modulated  signal  is  selected  as  64.0  times  the  modulating-tone 

' frequency,  the  results  from  the  FMDF00  program  including  table  5 are  valid 

for  any  values  of  the  modulated-signal  center  frequency. 

It  is  expected  that  no  even-order  distortion  will  be  produced  when  the 
bandpass  filter  is  tuned  to  the  modulated  signal.  Table  5 indicates  that  D2 
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and  D4  are  greater  than  260  dB  for  DFC  = 0.0.  This  result  can  be  interpreted 
to  mean  that  our  program  simulates  a laboratory  test  in  which  the  background 
noise  level  in  measurements  of  harmonic  distortions  is  at  least  260  dB  below 
the  desired  output  signal  level  (or  the  fundamental -frequency  component  of 
the  output  signal). 


5.4o  Third  Example  — FM  Noise 

As  the  third  and  last  example,  we  simulate  reception  of  an  FM  signal  in 
the  presence  of  white  Gaussian  noise.  This  is  a classic  problem.  As  in  the 
first  example  described  in  section  5.2,  we  simulate  reception  of  two  desired 
signals  with  and  without  modulation,  determine  two  corresponding  psophometri- 
cal ly-weighted  powers,  and  calculate  the  output  signal -to-noise  ratio  (SNR) 
as  the  ratio  of  the  difference  between  the  two  output  powers  to  the  output 
power  without  modulation.  Programming  the  core  of  the  program  for  this  ex- 
ample is  very  similar  to  programming  the  first  example  described  in  5.2.  The 
main  difference  between  this  example  and  the  first  example,  other  than  the 
differences  in  modulator,  signal  generator,  and  demodulator,  is  the  inclusion 
of  a bandpass  filter  in  this  example.  Since  the  interfering  signal  is  noise 
in  this  example,  the  result  largely  depends  on  the  RF  (or  IF)  bandwidth;  in 
other  words,  inclusion  of  a bandpass  filter  is  essential  in  this  simulation. 

The  FMWG00  program  that  is  written  for  illustration  purposes  is  listed 
in  table  6.  In  this  program,  we  have  selected 

NFC  = 128, 

UFS  = 500  (Hz), 

FC  = 30000  (Hz), 

where  NFC  is  the  number  of  frequency  components  in  the  spectrum  of  a signal, 
UFS  is  the  unit  frequency  spacing  in  the  spectrum,  and  FC  is  the  center  fre- 
quency of  the  modulated  signal.  We  have  also  assumed  a baseband  bandwidth  of 
3 kHz,  modulation  with  a single  tone  of  frequency  1 kHz,  and  a fourth-order 
bandpass  filter  with  the  Carson  bandwidth  and  with  no  ripple  in  the  passband. 
The  FMWG00  program  performs  simulation  for  nine  values  of  peak  frequency  de- 
viation and  for  11  values  of  the  input  signal -to-noise-density  ratio  (SNDR), 
which  is  defined  as  the  ratio  of  signal  power  (in  watts)  to  noise  power  den- 
sity (in  watts/hertz  or  joules).  Since  the  unit  frequency  spacing  (UFS)  is 
500  Hz  in  this  program,  this  program  simulates  measurements  of  output  SNR 
in  an  interval  of  2 ms  (or  2 cycles  of  the  modulating  signal).  This  short 
interval  of  measurement  will  result  in  fluctuations  in  the  output  SNR  values 
and  therefore  irregularities  in  the  curves  for  the  output  SNR  values  vs,  the 
input  SNR  (or  SNOR)  values. 

The  result  of  the  FMWG00  program  is  shown  in  table  7,  The  computation 
time  required  to  run  the  FMWG00  program  on  the  CDC-6600  computer  was  approxi- 
mately 15  seconds.  In  table  7,  FDV  denotes  the  peak  frequency  deviation  in 
Hz.  The  output  SNR  values  given  in  table  7 exhibit  some  irregularities,  es- 
pecially in  the  left  portion  of  the  table;  the  output  SNR  should  be  a mono- 
tonic increasing  function  of  the  input  SNDR  for  a fixed  FDV  value  theoreti- 
cally, while  the  output  SNR  shown  in  table  7 is  not  monotonic.  This 
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Table  6.  Fortran  Listing  of  the  FMWG00  Program 


PROGRAM  FMWGOe (OUTPUT, T A P E6=0UTPUT ) 

C FK  system  with  white  GAUSSIAN  NCISE 

C THE  INPUT  SIGNAL-TO-NOISE-OENSITY  RATIO  IS  THE  RATIO  OF  SIGTAL 
C POWER  (WATTS)  TO  NOISE  POWER  CENSITY  ( WA TTS/HERT Z) . 

C THE  OUTPUT  SIGNAL-TO-NOISE  RATIO  IS  CALCULATED  AS  THE,  RATIC  CF 
C OUTPUT  POWER  WITH  MODULATION  TO  OUTPUT  POWER  WITHOUT  VcDULA- 
C TION  MINUS  ONE. 

C DECLARATION  STATEMENTS 

DIMENSION  FC V(9)  ,SN0R(11)  ,SNRO(ll) 

COMPLEX  SFl  (la?)  ,SF«»  (129)  ,SF5  (128)  , 

1 VF1(12B),VF 2(128), VF3(128),VF4 (126), yF5(128) 

DATA  NFO/9/,  NSNR/11/ 

DATA  NFC/128/,  UFS/500.0/,  FC/30000.0/, 

1 FD\//100CO.,5000.,3000.,20  00.,1000.,5Q0..30Q.,20Q., 

2 100./, 

3 SFl/2* (( C.C, 0.0) ), (1.0,0. 0) ,125* ( (0. 0, 0.0)  ) / 

DATA  NOB/4/,  RPB/0.0/ 

DATA  SNRCMN/1 .02E-10/ 

DATA  NAME/6HFHWG00/,  LUN/6/ 

C DATA  PREPARATION 

10  DO  11  ISNR=1,NSNR 

SN0R(ISNR)=2*(ISNR-1)  +4  0 

11  continue 

C SIMULATION 

C - GENERATES  UMMOOULATEO  SIGNAL. 

20  CALL  MOOFH(NFC,UFS,SF1,FC,0.0,VF2) 

C - GENERATES  WHITE  GAUSSIAN  NOISE. 

CALL  SGWGN(NFr,UFS,2,NFC,\/F3) 

DO  89  IFO=l,NFO 
C - GENERATES  AN  FM  SIGNAL. 

FOI=FOY  (IFO) 

CALL  MOCFM( NFC,UFS,SF1,FC,FDI,VF1) 

C - ADJUSTS  THE  FILTER  BANDWIDTH  TO  THE  CARSCN  BANDWIDTH. 
3W3=2.0*(3000.0+F01) 

DO  59  ISNR=1,NSNR 

C - COMBINES  THE  FM  SIGNAL  WITH  THE  NOISE. 

CALL  C9CCT2(NFC, VFl ,VF3,SNOR(ISNR)  ,VF4) 

C - APPLIES  THE  COMPOSITE  SIGNAL  TO  THE  BANCPASS  FILTER  INFLT. 

CALL,  PPFCHtNFC  ,UFS  , NOB  , RPB,  FC  , BWB  , VF4,  OF 4) 

C - OEMOOULATES  THE  FILTER  OUTPUT  SIGNAL. 

CALL  CEMFM(NFC,UFS, VF4,FC,SF4) 

C - APPLIES  THE  OEMCDULATEO  SIGNAL  TO  THE  PSCPHOMETER  INPUT. 

CALL  PSPHCT (NFC,UFS ,SF4,PS4) 

C - COMBINES  THE  UNMODULATED  SIGNAL  WITH  THE  NOISE. 

CALL  C8CCT2,(NFC,  VF2  ,VF3,SNOR(  ISNR)  ,VF5) 

C - AOPLIES  THE  COMPCSITE  SIGNAL  TO  THE  BANDPASS  FILTER  INPUT. 

CALL  BPFCH (NFC.UFS, N0e,RPB,FC,BWB,VF5, VF5) 

C - OEMOOULATES  THE  FILTER  OUTPUT  SIGNAL, 

CALL  0EMFM(NFC,UFS, VF5,FC,SF5) 

C - APPLIES  the  demodulated  SIGNAL  TO  THE  PSOPHOHETER  INPUT. 

CALL  FSFHCT(NFC,UFS,SF5,PS5) 

C CALCULATION 

IF(PS5.LE.C.C)  GO  TO  42 
SNR0I=PS4/PS5-1. 0 

IF(SNROT.IE.SNROMN)  SNROI=SNROMN 
GO  TO  43 

42  SNROI=l.aE+20 
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Table  6.  Continued 


t 

I 

p 

I 

I 

* 
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Table  7.  Computer  Printout  Produced  by  the  FMWG00  Program, 
Listed  in  Table  6 

FM  SYSTEM  WITH  WHITE  GflUSEIAM  NOISE 

BASE3AN0  dANOWIOTH  = 302C  H? 

MODULATING  SIGNAL  A 1000-HZ  TONE 


OUTPUT  SIGnAL-TO-NOISE  RATIO  (36.  i 

INPUT  SNOR  <OB)  = I 

41.0  42.0  44.0  4f>.a  48.0  GO.O  52.0  54.0*  56.  0 58.  0 60.0  \ 

\ 

( 

7.4  8.0  5.6  1h.8  24.3  29.9  32.6  34.9  37.0  39.0  41.1  • 

8.6  16.0  17.9  20.3  22.7  24.9  27.0  29.1  31.1  33.2  35.2  | 

T.O  10.7  13.4  15.3  19.1  20.2  22.4  24.5  26.5  28.6  30.6  | 

j 

5.7  7.8  10.6  12.7  14.9  16.9  19.0  21.1  23.1  25.2  27.2  j 

-1.6  -.1  2.9  5.5  8.0  1C. 3 12.6  14.7  16.9  19.0  21.1 

5.7  -1.1  .1  2.0  4.0  6.0  8.0  9.9  11.9  13.8  15.8 

-4.7  -17.3  -99.9  -8.6  -3.5  -.2  2.5  4.8  7.0  9.1  11.2 

-13.6  -14.0  -13.1  -10.5  -7.5  -4.7  -2.1  .4  2.7  4.9  7.1 

17.8  -99.9  -99.9  -99.9  -99.9  -24.1  -9.3  -5.4  -2.6  -.3  1.9 


4 3 SNRO ( ISMR) =10. 0* ALOGIO  ISNROI) 

59  continue 

C PRINTING  OF  THE  RESULTS 

60  IF(IFO.EQ.l)  WRITE  (LUN, 20601  NAME, SNOR 

WRITE  (LUN, 2061)  FDI.SNRC 

89  CONTINUE 
CALL  EXIT 

C FORMAT  STATEMENTS 

2060  FORMAT(lHl,A6,ex,35HFM  SYSTEM  WITH  WHITE  GALSSIAN  NCISE// 

1 15X.28HBASEBAND  BANDWIDTH  = 3000  HZ/ 

2 15X, 38HMOOULATING  SIGNAL  A 1000-HZ  TONE//// 

3 27X.33HOUTFUT  S IG NAL -T C-NO ISE  RATIC  (C3)  / 

4 3X,3HFDV,3X,17HINPUT  SNDR  (03)  =/ 

5 iX,llF6.1/) 

2061  FORMAT ( IX /IX ,F6.0,2X,11F6.1) 

END 


irregularity  is  considered  to  stem  from  the  choice  of  a large  UFS  value,  as 
discussed  in  the  preceding  paragraph.  It  is  expected  that  this  irregularity 
can  be  reduced  by  taking  a small  UFS  value  (and  thus  increasing  the  measure- 
ment interval)  but  with  a corresponding  increase  in  the  computation  time.  Al- 
though the  center  frequency  of  the  modulated  signal  is  selected  as  30  kHz  in 
this  program,  the  result  shown  in  table  7 is  independent  of  the  modulated- 
signal  center  frequency. 


6.  CONCLUDING  REMARKS 

We  have  started  with  elementary  concepts  concerning  computer  simulation 
of  communication  systems  and  a simulation  model,  derived  guidelines  for  de- 
veloping a simulation  model,  described  a model  that  has  been  developed  along 
these  guidelines,  and  explained  how  to  use  the  model.  The  model  described  in 
this  report  is  still  under  development.  For  further  development  of  the  model, 
we  make  the  following  remarks: 

The  model  presented  here  includes  modulators  and  demodulators  of  the  AM, 
FM,  FSK,  phase-modulation  (analog),  and  SSB  systems.  Because  of  its  increased 
use,  inclusion  of  a phase-shift-keying  (PSK)  system  is  highly  desirable.  It 
is  anticipated  that  inclusion  of  PSK  will  necessitate  simulation  of  phase- 
lock  loop  (PLL),  which  has  a wide  field  of  application. 

At  the  present  time,  no  baseband  signal  generator  is  included  in  the 
model.  The  addition  of  generators  of  some  typical  baseband  signals,  such  as 
a voice  signal,  will  be  desirable. 

Only  a Gaussian  noise  generator  is  included  in  the  model  as  a noise  gen- 
erator. Inclusion  of  non-Gaussian  noise  generators  for  atmospheric  noise  and 
man-made  noise  is  expected  to  make  the  model  more  versatile.  It  is  also  con- 
sidered desirable  to  include  interference  generators  that  simulate  uninten- 
tional emissions  from  industry  or  various  types  of  home  equipment. 

Only  a Chebyshev-type  bandpass  filter  is  included  in  the  model.  It  is 
expected  that  inclusion  of  an  elliptic-function-type  filter  will  make  the 
model  more  versatile  and  better  able  to  simulate  the  transfer  characteristics 
of  particular  receiving  systems. 

In  addition,  techniques  for  specifying  more  specific  and  detailed  per- 
formance measures  (such  as  a’^ticulation  scores  for  voice  signals,  bit  error 
rates  for  digital  signals,  direction  finding  errors,  etc.)  need  to  be  devel- 
oped. 

Inclusions  of  new  types  of  modulators  and  demodulators,  signal  and  noise 
generators,  bandpass  and  lowpass  filters,  etc.,  described  above  can  be  made 
by’adding  to  the  model  additional  subroutines  that  simulates  these  system  com- 
ponents. Those  subroutines  already  programmed  and  presented  in  this  report 
need  no  modifications. 

Articulation  index  (AI)  is  sometimes  used  as  a measure  of  voice  intelli- 
gibility (ANSI,  1970).  In  our  model,  however,  we  have  intentionally  omitted 
AI  because  the  validity  of  AI  for  radio-interference  problems  is  questionable. 
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APPENDIX  A 

FORTRAN  LISTINGS  OF  THE  MODEL 


Complete  Fortran  listings  of  all  subprograms  included  in  the  model  are 
given  in  this  appendix  in  alphabetical  order„  All  subroutines  are  written  in 
ANSI  Standard  Fortran o 

The  RANF  function  called  by  some  subroutines  is  a random-number  genera- 
tor, Repeated  calls  to  RANF  generates  a sequence  of  random  numbers  uniformly 
distributed  in  a range  extending  from  0 to  I,  If  the  computer  system- library 
does  not  include  this  function,  the  user  of  the  model  must  supply  it. 

In  addition,  the  following  intrinsic  and  basic  library  functions  are 
called  by  some  subroutines.  They  are  included  in  most  computer  systems  that 
accept  ANSI  Standard  Fortran, 


Description 


CMPLX 

FLOAT 


ATAN2 


Absolute  value  of  a real  number. 

Integer  part  of  a real  number  (truncation) 
Real  modulo  function. 

A complex  number  from  two  real  numbers. 
Conversion  from  integer  to  real. 

Natural  logarithmic  function. 

Arctangent  function. 

Cosine  function. 

Exponential  function. 

Sine  function. 

Square-root  function. 


SUtiROUTINc  BPFCM  tNFC,UFS»NOB.KPa,FCB,BHB,  VFlf  \/F2) 

C This  subroutine  simulates  a bandpass  filter  (BPF)  of  chebyshev 

C TYPE. 

C THE  INPUT  PARAMETEsS  ARE 

C NFC  = NUMBER  OF  FREQUENCY  COMPONENTS 

C IMUST  BE  POSITIi/E  BUT  NOT  CREATER  THAN  2O<«0)i 

C UFS  = UNIT  FREQUENCY  SPACING  (MUST  BE  POSITIVE), 

C . NOB  = ORCER  CF  THE  OFF 

C (MUST  NOT  oE  NEGATIVE  OR  GREATER  THAN  1001, 

C RPB  = RIPPLE  (Ob!  RITHIN  THE  FASS3AN0  OF  THE  BPF 

C (MUST  BE  NONNEGATIVt! , 

C FCB  = CENTER  FREQUENCY  OF  THE  BPF 

C (MUST  BE  POSITIVE), 

C BWB  = 3-DB  BANDWIJIH  OF  THE  BPF  (MUST  BE  POSITIVE), 

C VFl  = COMPLEX  ARRAY  CF  UIMENSION  NFC  CONTAINING  THE 

C SPECTRUM  COMPONENTS  OF  THE  INPUT  VOLTAGE. 

c THE  OUTPUT  Parameter  is 

C VF2  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE 

C SPECTRUM  COMPONENTS  OF  THE  OUTPUT  VOLTAGE 

C ARE  TU  3£  SIOREO. 

C THIS  subroutine  CALLS  THE  FRBPFC  SUBROUTINE, 
c declaration  statements 

complex  VFl (NFC) , VF2 (NFC) 

COMPLEX  MF9(2OM0) 

data  NFCMX/20L!)/,  NOBMX/lOO/ 

OATA  NFCPV/0/,  UFSPV/0. j/,NOBPV/-l/,  RP8PV/-1.0/, 

1 FCBPV/0.:/,  BHBPV/O.u/ 

OATA  LUN/b/ 

C PRELIMINARY  PROCESSING 
IG  NFC0=NFC 
UFS0=UFS 
N080=NOB 
RPB0=RP3 
FCB0=FCB 
BWBa=3WB 

IF(NFCa.LE.C.OR.NFCO.GT.NFCMX)  GO  TO  90 

IF(UFSO.LE.O.O)  GO  TO  90 

IFINOBO.LT.O.OR.NJBQ.GT.NOGMX)  GO  TO  90 

IF(RPB3.LT.0.I)  GO  TO  90 

IF(FC  .LT.0.0)  GO  TO  90 

IF(8WBQ.LT.0.0)  GO  TO  90 

IF(N030.EQ.O)  GO  TO  Ru 

C CALCULATES  FREQUENCY  RESPONSE  OF  THE  BPF  WHEN  NECESSARY. 

20  IF  (NFCO.EQ.NFCPV.AND.UFSC.EQ.UFSPV.ANO. 

1 NO30.EQ.NO3FV.ANO.RPB0.CQ.RPBPV. ANO. 

2 FC30.EQ.FC3PV.ANO.BHBJ.EQ.3WBPV)  GO  TO  30 

CALL  FRBPFC(NOB0,RP8u,FC30,3W80,UFSu,l,NFC0,HF8) 

C CALCULATES  OUTPUT  SPECTRUM. 

C - N08.NE.O 

30  DO  31  I = l,NFCCi 

VF2  (I ) =HF8(I) *VF1 (I ) 

31  CONTINUE 
RETURN 

C - NOB.EQ.Q 

40  00  41  I^l.NFCb 

VF2(I)=VF1II) 

41  CONTINUE 
RETURN 

C ERROR  EXIT 

90  WRITE  (LUN,209C)  NFCC ,UFS 0 , NO80 , RP30 , FCBO , BWBO 
RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 
2090  F0RMAT(1X/31H  *•*  IMPROPER  INPUT  VALUE(S)./ 

1 3X.5HNFC  =, Ib,5X,5HUFS  = , E 11 . 3 , 5 X , 5HN0B  =,I6, 


GO  TO  90 
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£ 5X,i>HRPB  =f€ll.  3,5X,5HFC6  =,E11.  3,5X  ,5H8W8  *,tll.3/ 
3 35H  ERROR  OtTECTEO  IN  ROUTINE  BPFCH  /) 

END 


SUBROUTINE  C8CCT2  ( NFC  t VFl  t \/F2,R0B,  VP 31 
THIS  SUBROUTINE  SIMULATES  A RAO lO-FkEUUENCV  IRF)  COMBINING 
CIRCUIT  THAT  linearly  COMBINES  TWO  RF  SIGNALS  WITH  A SPECI 
FIcD  RATIO. 

THE  LEVEL  OF  THE  FIRST  INPUT  SIGNAL  IS  UNCHANGED  IN  THIS 
CIRCUIT. 

THE  INPUT  PARAMETERS  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS  IN  THE  SPECTRUM 
OF  EACH  INPUT  SIGNAL. 

VFl  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 
VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  FIRST  INPUT 
SIGNAL. 

VF2  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 

VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  SECOND  INPUT 
SIGNAL. 

ROB  = SIGNAL-TO-INTERFERENCE  RATIO  OR  SI GNAL-TO-NOISE- 
OENSITY  RATIO  BETWEEN  THE  FIRST  AND  THE  SECOND 
INPl/T  SIGNALS  IN  08. 

THE  OUTPUT  PARAMETER  IS 

VF3  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  COMPOSITE  SIGNAL  ARE 
TO  BE  STORED. 

DECLARATION  STATEMENTS 

DIMENSION  VFl (EmJ .VF2J64) . VF3<64» 

DATA  RPV/0.0/.  VR/1.3/ 

PRELIMINARY  PROCESSING 
10  NO=NFC 
RO  = RDB 
NT2=N0*2 

CONVERTS  SIGNAL-TO-INTERFERENCE  RATIO  OR  SIG NAL-TO-NOISE- 
OENSITY  RATIO  FROM  DB  TO  VOLTAGE  RATIO  WHEN  NECESSARY. 

20  IF(RO.EQ.RPV)  GO  TO  30 
RPV=RO 

VR=EXP< -0.11512 925M6»R0> 

LINEARLY  COMBINES  THE  SPECTRA  OF  THE  TWO  SIGNALS  ON 
COMPONENT-BY-COMPONENT  BASIS. 

30  DO  31  D=1.NT2 

VF3IJ)=VF1(D) ♦VR*VF2ID) 

31  CONTINUE 
RETURN 
END 


SUBROUTINE  CBCCT3 (NFC. VFl . VF2 . VF 3,R1 20 3 .R13DB. VF4) 

THIS  SUBROUTINE  SIMULATES  A RADIO-FREQUENCY  IRF)  COMBINING 
CIRCUIT  THAT  LINEARLY  COMBINES  THREE  RF  SIGNALS  WITH  SPECI- 
FIED RATIOS. 

The  level  of  the  first  input  signal  is  unchanged  in  this 

CIRCUIT. 

THE  INPUT  PARAMETERS  ARE 

NFC  * NUMBER  OF  FREQUENCY  COMPONENTS  IN  THE  SPECTRUM 
OF  EACH  INPUT  SIGNAL. 

VFl  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 
VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  FIRST  INPUT 
SIGNAL. 

VF2  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 

VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  SECOND  INPUT 
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I C SIGNAL, 

f C VF3  = complex  Ar<KAY  OF  DIMENSION  NFC  CONTAINING  THE 


f C \/OLTAGt  SPECTRLI"^  COMPONENTS  OF  THE  THIRO  INPUT 

C SIGNAL, 

C R12DB  = SIGNAL-TO-INTERFERENCE  RATIO  OR  SIGNAL-TO- 

C NOISE-OENSITY  RATIO  3ETWEEN  THE  FIRST  AND  THE 

C SECOND  INPUT  SIGNALS  IN  OB, 

C'  R13QB  - SIGNAL-TO-INTERFERlWCE  RATIO  OR  SIGNAL-TO- 
■ C NOISE-CENSITY  RATIO  BETWEEN  THE  FIRST  AND  THE 

c Third  input  signals  in  do. 

c the  output  parameter  is 

c VF4  = COMPLEX  ARRAY  OF  OIMcNSION  NFC  WHERE  THE 

C VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  COMPOSITE 

C SIGNAL  ARE  TO  BE  STORED. 

C DECLARATION  iTATEMENTS 

, DIMENSION  VFl 164) , VF2t 64) , VF3 (b4) , VF4 < 64) 

DATA  R2PV/0.0/,  R3PV/0.0/,  VR2/1.0/,  VR3/l.a/ 

C PRELIMINARY  PROCESSING 
/ 10  NC=NFC 

R2=R120B 

R3=R13JB 

NT2=N0*2 

C CONVERTS  SIGNAL-TO-iNTERFERENCE  RATIOS  AND/OR  SI GNAL-TO-NOISE- 
C DENSITY  -RATIOS  FROM  03  TO  VOLTAGE  RATIOS  WHEN  NECESSARY. 

20  IF (R2.EQ.R2PV)  GO  TO  21 

R2PV-R2 

I VR2=EXP1-0.1131232546»R2) 

21  IF1R3.EU.R3PV)  GO  TO  30 

R3PV=R3 

VR3=EXP 1-0. 115129254b*R3) 

C linearly  combines  the  SPECTRA  OF  THE  THREE  SIGNALS  ON 
C COMPONENT-3Y-GOMPONENT  BASIS. 

30  DO  31  J-1,NT2 

VF4(J)=VF1(J)4VR2*VF2(J)+VR3*VF3(U) 

31  CONTINUE 
RETURN 
END 


SUBROUTINE  DEMAM  INFC,VF,SF) 

C THIS  SUBROUTINE  SIMULATES  AN  AMPLI lUDE-MOOUL AT  I ON  (AM)  DEMODU- 
C LATOR. 

C THE  INPUT  PARAMETERS  ARE 

C NFC  = NUMBER  OF  FREQUENCY  COMPONENTS. 

C VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 

C VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  DEMODULATOR 

C INPUT  SIGNAL. 

C THE  OUTPUT  PARAMETER  IS 

C SF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 

C SPECTRUM  COMPONENTS  OF  THE  DEMODULATED  SIGNAL 

C ARE  TO  BE  STORED. 

C THIS  SUBROUTINE  CALLS  THE  OFT  SUBROUTINE. 

C DECLARATION  STATEMENT 

DIMENSION  VF164) ,SFI64) 

C PRELIMINARY  PROCESSING 
10  N0=NFC 
NT2=N0*2 

C FOURIER-TRANSFORMS  THE  INPUT  SIGNAL  FROM  FREQUENCY  DOMAIN  TO 

c time  domain. 

20  CALL  OFT(NO,VF,1,SF) 

C EXTRACTS  THE  AMPLITUDE  INFORMATION  AT  EACH  SAMPLING  POINT  IN 
C the  time  DOMAIN  AND  STORES  THE  RESULTS  IN  THE  SF  ARRAY. 

30  00  31  J=2.NT2,£ 
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SF< J-1)=SQRT{SF( J-1)**2»SF (J) **2) 

SF(J)  =0.0 

31  continue 

C FCURIER-TRANSFORMS  THE  OENOOULATEO  SIGNAL  FROM  TIME  DOMAIN  TO 
C FREQUENCY  DOMAIN, 

AO  CALL  OFT<N0»SF,-l,SF) 

C ADJUSTS  THE  OEMOOULATEO'SIGNAL  SPECTRUM  IN  SUCH  A HAY  THAT 
C-ITS  UPPER  HALF  BECOMES  ZERO.  , 

50  J2=NT2*2 

00  51  J1=A«N0.2 

J2=J2-2 

SF  ( Jl-1)  =SF  ( Jl-U  *SF  U2-1) 

SF(Jl)  =SF(J1)  -SF«J2> 

SF<J2-1J  =0,0 
SF(J2)  =0.0 

51  continue 
RETURN 
END 


SUBROUTINE  OEMFM  « NFC, UFS , VF ,FCO, SF) 

THIS  SUBROUTINE  SIMULATES  A FREQUENCY-MODULATION  (FM)  OEMOOU- 
LATOR. 

THE  INPUT  PARAMETERS  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS, 

UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 

VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 
voltage  spectrum  COMPONENTS  OF  THE  DEMODULATOR 
INPUT  SIGNAL, 

FCO  = CENTER  FREQUENCY  OF  THE  OEMOOULATOR. 

THE  OUTPUT  PARAMETER  IS 

SF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  OEMOOULATEO  SIGNAL 
ARE  TO  BE  STORED. 

THIS  SUBROUTINE  CALLS  THE  OFT  SUBROUTINE. 

DECLARATION  STATEMENTS 

DIMENSION  VF(6<,)  ,SF(6<»» 

C PRELIMINARY  PROCESSING 
10  NC=NFC 
U0=UFS 

F0=U0* AINU FCO/ U0+0.5> 

NT2=N0*2 

C DIFFERENTIATES  THE  INPUT  SIGNAL  IN  THE  FREQUENCY  DOMAIN  AND 
C STORES  THE  RESULT  IN  THE  SF  ARRAY. 

20  FI=-U0 

DO  21  J=2,NT2,2 

FI=FI+U0 

SF(J-1)*-FI*VFCJ) 

SF(J)  = fI*VF<J-U 

21  CONTINUE 

C FOURIER-TRANSFORMS  THE  SF  AND  VF  ARRAYS  FROM  FREQUENCY  DOMAIN 
C TO  TIME  DOMAIN. 

30  CALL  OFT(NO,SF,1,SF) 

CALL  0FTtN0,VF,l,VF) 

C CALCULATES  THE  INSTANTANEOUS  FREQUENCY  AT  EACH  SAMPLING  POINT 
C IN  THE  TIME  DOMAIN  AND  STORES  THE  RESULTS  IN  THE  SF  ARRAY. 

40  DO  41  J«2,NT2,2 

A=SF<J-1» 

9«SFU) 

C»VF(J-1) 

0=VF( J) 

SF(J-l) CB*C-A*D)/«C*C»D*D)-FO) 

SF(J)  =0.0 
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41  CONTINUE 

C FOURIER-TRaNSFOkMS  THE  OtHOOULaTEO  SIGNai  FSCM  TIME  OOM£IN  TO 
C FRtCUtNCY  OOMaiN. 

50  CaiL  OFT (N0,SF,-1, SF) 

C RESTORES  THE  INPUT  \lf  ARPar. 

CALL  OFT  (NO^VF,-!,  \/F) 

C adjusts  the  GtMCuULaTEJ-SIGNaL  SPECTRUM  IN  SUCH  a WAY  Th^T 
CITS  UPPER  HALF  BECOMES  ZERO. 

60  J2=NT2tZ 

DO  61  Jl=4,N0t2 

J2=J2*2 

SF  ( Jl-1)  = SF  ( Jl-l ) 4-SF  ( J2-1) 

SFIJll  =SF(J1)  -SF«J2I 

SF(J2-l)=C.O 

SF(J2)  =C.O 

61  CONTINUE 
RETURN 
ENO 


SUBROUTINE  OEMFS<(NFC,UFS,VF,Fco,NI3,IB) 

THIS  SUBROUTINE  SIMULATES  A F Rt QUENC Y -SH I FT -KE Y I NG  (FSKl  DE- 
MODULATOR. 

THE  INPUT  PARAMETERS  ARE 

NFC  = NUMBER  OF  FREOUENCY  COMPONENTS 
IMUST  NOT  EXCEtO  2048)  . 

UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 
y/F  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 
VOLTAGE  spectrum  COMPONENTS  OF  THE  DEMODULATOR 
INPUT  SIGNAL, 

FCD  = CENTER  FREQUENCY  OF  THE  DEMODULATOR, 

NIB  = NUMBER  OF  INFORMATION  BITS  IN  A PERIOD. 

THE  OUTPUT  PARAMETER  IS 

10  = ARRAY  OF  DIMENSION  NIB  HHERE  THE  DEMODULATED 

INFORMATION  BITS  ARE  TO  BE  STORED. 

THIS  SUBROUTINE  CALLS  THE  OFT  SUBROUTINE. 

DECLARATION  STATEMENTS 

DIMENSION  VF(6hI  ,I‘B(lb) 

DIMENSION  HK(-*oB6) 

DATA  NFCMX/2348/ 

DATA  LUN/6/ 

C PRELIMINARY  PROCESSING 
10  N0=NFC 

IFINO.GT.NFCMX)  GO  TO  90 

U0=UFS 

F0=U0*AINT (FCD/U0M0.5) 

NlsNIB 

NT2=N0*2 

NI2=NI*2 

C OIFFERENTIATES  THE  INPU»  SIGNAL  IN  ImE  FREQUENCY  DOMAIN  AND 
C STORES  THE  RESULT  IN  THE  WK  ARRAY. 

20  FI=-U0 

00  21  J=2,NT2,2 

FI=F1*U0 

MKU-1)  = -FI*VF(J) 

WK(J)  = FI»WF(J-1I 

21  CONTINUE 

c fourier-transforms  the  vf  and  wk  arrays  from  frequency  domain 
C TO  time  domain. 

30  CALL  OFTINC,VF,l,VF) 

CALL  OFUNO,WK,1,HK) 

C calculates  the  instantaneous  FREuUEnOY  AT  each  sampling  point 
c IN  the  time  domain  and  stores  the  results  in  The  hk  array. 
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40  CO  41  1=1, NO 

J=I»I 
A=HK« J-1) 

B=WKIJ) 

C=VF(J-1) 

0=VF(J) 

HK<I)=(a*C-A»O)/(C*C»O*O»-F0 
• 41  CONTINUE  ' 

C RESTORES  THE  INPUT  VF  ARRAY. 

50  CALL  OFT(Nu»VF,-l,YF) 

C SIMULATES  3IT  SYNCHRONIZATION. 
bu  NS=N0/NI2+1 
SHMX=0.0 
00  64  IS=1,NS 
J=NC^IS 
NK(Jt=MK(ISi 
SM=0.0 
IIBPY=1 
J = I5 

S3=MK(J) 

00  63  1=2, NC 

II3=t (I-1»*NI) /N0*l 
IF(IIB.EQ.IIBPY)  CO  TO  62 

61  SM=SH«ABS (SB) 

IIBPV=IIB 

SB=0.C 

62  J=J*1 
S3=SB*HKIJ) 

63  CONTINUE 
SM=SM4ABS«SB> 

IF(SM.LE.SMMX)  GO  TO  64 

SMMX=SM 

ISMX=IS 

64  CONTINUE 

C estimates  the  transmitted  bit  SEQUENCE  AND  STORES  THE  RESULTS 
C IN  THE  IB  ARRAY. 

65  IIBPY=1 
J=ISMX 
SB=HK(J> 

DO  6S  1=2. NO 

IIB=(  (I-l)  •ND/NO  + l 
IFlIIB.EQ.IIBPVi  GO  TO  67 

66  I3(IIBPV)=G 

IF (SB. GT. 0.3)  IB(II3PV)=1 

IIBPY=IIB 

SB=0.0 

67  J*J*1 
S3=SB»HK(J) 

6S  CONTINUE 
IB(II3)=0 

IFCSB.GT.C.C)  IBdIBXl 

return 

C ERROR  EXIT 

90  WRITE  (LUN,2090)  NO 
RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 
2090  FORMATI1X/31H  IMPROPER  INPUT  VALUE  IS)./ 

1 6H  NFC  =.16/ 

2 35H  ERROR  DETECTED  IN  ROUTINE  OEMFSK/I 
END 
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SUBRCUTINt  btMFh.'l  (Nf  C,UFS.  VFtSF) 

THIS  SU3R0JTINt  SIMULftTES  A PHASt-MOOULAT ION  DEMODULATOR. 
the  input  parameters  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS, 

UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 

VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 
YOLTAGE  SPECTRUM  COMPONENTS  OF  THE  DEMODULATOR. 
INPUT  SIGNAL. 

OUTPUT  PARAMETER  IS 
SF  = COMPLEX  ARRAY  CF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  DEMODULATED  SIGNAL 

THE  OUTPUT 
WILL  EQUAL 

A SINGLE  SINUSOID  WITH 
THE  CF  CONSTANT  IN  THE 
AND  THE  CH  IS  ONE  hAlF 
THIS  SUBROUTINE  CALLS 
DECLARATION  STATEMENTS 
DIMENSION  VF(E«.) 

DATA  TWOPI/6.2831 
1 CH/4. ALEeSOM 

PREL 
13 
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C 
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c 
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8857bS8 

381 

5837/ 

IMINARY  PROCESSING 
N3=NFC 


U0=UFS 

NT2=N0*2 

FLN2=NT2 

C FOURIER-TRANSFORMS  THE  INPUT  SIGNAL  FROM  FREQUENCY  DOMAIN  TO 
C TIME  DOMAIN. 

33  CALL  DFT(N0,VF,1,SF» 

C EXTRACTS  THE  PHASE  INFORMATION  AT  EACH  SAMPLING  POINT  IN 
C THE  TIME  DOMAIN  AND  STORES  THE  RESULTS  IN  THE  SF  ARRAY. 

AO  DO  A2  J=2,NT2,2 
CYCL=0.0 


A=SF( J-1) 

3=SF ( J) 

IF<A.NE.O.C.OR.3.NE.O.OI  CYCL=ATAN2(B,A)/TW0PI 

IF(J.EQ.2>  GO  TO  A1 

DCYCL=CYCL-CYCLPV 

IF (DCYCL.LT .C .0)  CYCL=CYCL+AINT(-DCYCL*1.0) 

Al  SF<J-1)=CYCL 
SF(J)  -0.0 
CYCLP  =CYCL 
A2  CONTINUE 

AS  DCYCL  = AINTJSF(NT2-l)-SFtlltl.(l) 

0CYCL=DCYCL/FLN2 
DO  Ab  J=2,NT2,2 

SF  (J-H=CF*  (SF  (J-ll  -OCYCL*FLOAT  (J-2)  J 
Ab  CONTINUE 

C FOURIER-TRANSFORMS  THE  DEMODULATED  SIGNAL  FROM  TIME  DOMAIN  TO 
C FREQUENCY  DOMAIN. 

so  CALL  OFT (N0,SF,-l,SF) 

C ADJUSTS  THE  DEMOOULATEO-SIGNAL  SPECTRUM  IN  SUCH  A HAY  THAT 
C ITS  UPPER  HALF  BECOMES  ZERO. 
bO  J2=NT2+2 

DO  bl  J1=A,NQ,2 


J2=  J2-2 

SF  (Jl-1)  =SF  ( Jl-l)  ♦SFU2-1) 
SF(Jl)  sSF(Jl)  -SF(J2) 
SFt  J2-U=0.0 
SF(J2I  =0.0 
bl  CONTINUE 

70  IF(SF(ll .GE.-CH)  GO  TO  7l 
SF(1)=SF(1)*CF 
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RETURN 

71  IF(SFU)  .LE.CHJ  RETURN 

SF(1)=SF(1)-CF 
RETURN 
END 


SUdROUTINE  0EMSS3(NFC,UFSt  \/F,FRutSF» 

This  subroutine  simulates  a SINGLE-SIOEBANO  (SSBI  oemodulator. 

THE  INPUT  parameters  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS* 

UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM. 

VF  = COMPLEX  ARkAY  OF  DIMENSION  NFC  CONTAINING  THE 

(/oltage  spectrum  components  of  the  oemodulator 

INPUT  SIGNAL* 

FRO  = REFERENCE  FREQUENCY  OF  THE  DEMOOULATOR. 

THE  OUTPUT  PARAMETER  IS 

SF  s COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  OEMOOULATEO  SIGNAL 
ARE  TO  BE  STORED. 

THIS  SUBROUTINE  CALLS  THE  RANF  FUNCTION. 

DECLARATION  STATEMENTS  j 

DIMENSION  VF  (bx)  *SF  (61* ) i 

DATA  TH0PI/6.E83185337179S9/  j 

C PRELIMINARY  PROCESSING  ! 

to  N0=NFC  j 

UG=UFS  j 

F0=FRD  I 

NT2=2*N0  j 

IFR=F0/U0*1.6  I 

IFRT2=1FR*2  | 

C GENERATES  THE  CARRIER  LOCALLY.  3 

20  PHI=THOPI* (RANF(0» -0.5»  | 

Cl=COS(PHI» 

C2=SIN(PHI)  i 

C TRANSLATES  EACH  SPECTRAL  COMPONENT  OF  THE  INPUT  SIGNAL  DORN  * 

C TO  THE  BASEBAND  AND  STORES  THE  RESULTS  IN  THE  SF  ARRAY.  ; 

30  J1=IFRT2»2  : 

J2=IFRT2-2  [ 

DO  3<*  K=2*NT2*2 

S1=0.0 
S2=0.0 
Jl=Jl-2 
J2=J2>2 


IFtJl.GT.NT2) 

GO 

TO 

33 

IF (K.EQ.2) 

GO 

TO 

32 

IF< Jl.LE.b) 

GO 

TO 

31 

V1=VF(J1-1) 

V2=VF(J1» 

Sl=  Cl*Vl+C2*V2 
S2=-C1*V2+C2*V1 

31  IFU2.GT.NT2)  GO  TO  33 

32  V1=VF(J2-1) 

V2=VF(J2) 

S1=S1*C1*W1»C2*V2 

S2»S2tCl*V2-C2*Vl 

33  SF(K-1)=S1 
SFOO  *S2 

3n  CONTINUE 
RETURN 
ENO 
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SUBROUTINE  0FT(N.X1,1S0N|X^> 

this  subroutine  performs  discrete  FOURIER  TRANSFORM.  THIS  IS 
A MODIFIEO  VERSION  FROM  THE  DISCRETE  COMPLEX  FAST  FOURIER 
TRANSFORM  PACKAGE  FROGRAMMEO  BY  L.  OAVIO  LEWIS  ANO  MARIE  WEST 
OF  THE  U.S.  OEPARTHENT  OF  COMMERCE  BOULDER  LABORATORIES.  tCF. 
SINGLETON,  COMMON.  ACM,  VOL.  U,  PP . 776-779,  1968,  AND  VOL. 
12,  P.  187,  1969) 

the  input  parametefs  are  f 

N = OlhENSION  OF  THE  XI  ANU  X2  ARRAYS 
= 2»»M  (M=U,1,2,...,1G), 

XI  = COMPLEX  ARRAY  OF  DIMENSION  N CONTAINING  THE 
INPUT  SEQUENCE  TO  DE  TRANSFORMED, 

ISGH  = SIGN  IN  THE  EXPONENT 

= ♦!  FROM  FREQUENCY  SPECTRUM  TO  TIME  SEQUENCE 
= -1  FROM  TIME  SEQUENCE  TO  FREQUENCY  SPECTRUM. 

THE  OUTPUT  PARAMETER  IS 


X2  = COMPLEX  ARRAY  OF  DIMENSION  N WHERE  THE 
TRANSFORMED  SEQUENCE  IS  TO  BE  STORED. 

The  P2  array  contains  powers  of  two  --  P2(I)=2**(I-1) 

THE  ST  ARRAY  IS  SINE  ARRAY  --  ST ( I ) =SIN (PI / ( 2** I) ) . 

declaration  statements 

DIMENSION  Xl(2j<,8)  ,X2(20>,8) 

DIMENSION  P2(15),STU5) 

INTEGER  P2 


REAL  IM 

DATA  P2(l) ,P2 (21 ,P2 (3) ,P2 (A) ,P2 (5) ,P2(6) ,P2 (7) ,P2t 8) , 
P2 (9) ,P2 1101 ,P2( 111 ,P2 (12) ,P2( 13) ,P2( 14) ,P2 (15)/ 
1, 2, M, 8, 16, 32, 6h, 128, 256, 512, 10 24, 2048, 40 96, 8192, 
16384/ 

DATA  ST(1),ST(2),ST(3),ST(4),ST(5),ST(6),ST(7),ST(8), 
ST(9)  ,ST  (10)  ,ST (11) ,ST (12) ,ST(13) ,ST(14) ,ST(15)/ 
.lOOOOQOOOOOOOOOE^l,  .70  710678118654  8E«'0, 
.3826834 32 3650 90E»0,  . 195 090 3220 1 612 8E«0 , 

•490676743274180E-1, 

. 122715382857199E-1, 
.306795676296598E-2, 
.7669903l8742705E-3i 
. 191747597310703E-3, 


1 

2 

3 

[ 

1 

2 

3 

4 

5 

6 
7 

7 

8 


DATA 
PRELIMINARY 
10  NQsN 


. 980171403295606E-1, 
.245412285229123E-1, 
.613588464915448E-2, 
. 15339801do28477E-2, 
.383495187571396E-3, 
.9587379903597732-4/ 
LUN/6/ 

PROCESSING 


DO  11  1=1,15 

.TF(N0.EQ.P2(I)  ) GO  TO  12 
11  CONTINUE 
GO  TO  90 


12  M0=I-1 
NT2=N0*t 
DO  13  K=1,NT2 


X2(K) =Xl (K) 

13  CONTINUE 

IF(MO.EQ.O)  RETURN 

NP1=N0«-1 

NM1=N0-1 

NQ=N0/4 

ISGN0=ISGN 

C REVERSIBLE  PERMUTATION  OF  THE  X2  ARRAY  FROM  NORMAL  SEQUENCE 
C TO  REVERSE  BINARY  SEQUENCE. 

20  IF(MG.EQ.l)  GO  TO  30 

I=M0 
K0  = 1 


KI  = 2 
KJsNMl 

21  K0=K0*P2(1) 
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GO  TO  22 


i 

i 


P2a>=-P2«i) 

IF1P2(1).GE.0>  GO  TO  23 

IF (Kl. bT.KO. OR.KJ.lt. Ka) 

KIT2=KI*2 
KuT2=K3*2 
RE  = X2  (KIT2-1) 

IM=X2 (KIT2I 
X2(KIT2-U=X2(K0T2-1I 
X2IKIT2I  =X2(KCT2) 

X2(KuT2-l)=RE 
X2(K0T2)  =IM 
IF(KJ.EQ.KO)  GO  TO  22 
KK=NP1-KJ 
KKT2=K<*2 
KJT2=KJ*2 
R£  = X2  (KKT2-1) 

IM=X2 (KKT2) 

X2IKKT2-1)=X2(KJT2-1) 

X2(KKT2»  =X2(KJT2I 
X2 (KJT2-1) =RE 
X2IKJT2)  =IM 

22  I=MO 
KI=K1*1 
KJ=K.:-1 
GO  1 ; 21 

23  IF(I.Ea.2)  GO  TO  30 

1 = 1-1 
GO  TC  21 

C FOURIER  TRANSFORM  FROM  INPUT  IN  REi/cRSE  BINARY  SEQUENCE  TO 
C OUTPUT  IN  NORMAL  SEQUENCE. 

30  00  31  KK=1,N0.2 

KKT2=KK*2 

RE=X2IKKT2-1)-X2IKKT2+1I 
IM=X2(KKT2J  -X2CKKT2*2I 
X2  <KKT2-1)=X2  <KKT2-1) ♦X2(KKT2»1) 

X2(KKT2»  =X2(KKT2»  ♦X^^K^T^♦2) 

X2 {KKT2*l)=Rfc 
X2 (KKT2+2»=IM 

31  CONTINUE 

40  IF(MO.£Q.l)  GO  TO  50 
EXPS=1.0 

IF(ISGNO.LT.C)  EXPSs-1.0 

JSPAN=1 
K=NQ 

DO  43  I=2.M0 

SO=-ST«I-l) 

CD  = 2.0*ST(l)*STa) 

R=-2. 0*CD 
CN=1.0 
CM=0.0 
SN=0.0 
SM=EXPS 
JSPANHsJSPAN 
JSPAN  sJSPANtJSPAN 
jj  = 0 
KK  = 1 

41  KS'KK+JSPAN 
KKT2=KK*2 

KST2=KS*2  ^ 

RE*CN*X2 IKST2-1) -SN*X2 JKST2) 

IM=SN*X2(KST2-1)  ♦CN*X2(ICST2» 

X2(KST2-1)=X2(KKT2-1)-RE 
X2IKST2>  *X2(KKT2)  -IM 
X2JKKT2-1)  = X2  (KKT2-H  *RE 
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X2(KKT2)  =X2(K<T2)  ♦!« 

KK=KK*JSPANH 

KS=KS*JSPANH 

KKT2=KK*2 

KST2=KS*2 

Rt  = CM*X2  t)CST2-l)  -SM»X2  (KST2) 

IM=SM*X2<KST2-U +CM»X2(KST2» 

X2  1KST2-U=X2  (KKT2-1I -RE  ' 

X2(KST2)  =X2tKKT2)  -IM 
X2  (KKT 2- 1 ) = X2  (K«;T2-1)  >RE 
X2(KKT2)  =X2(K<T2)  +IM 
KK=<S*JSPANH 

IFIKK.LT.NCi)  GO  TO  41 

KKsKK-NMl 

JJ=JJ4K 

IFCJJ.GE.MQ)  GO  TO  42 

CD=R*CN.CD 

CN=CO*CN 

SM=CN*£XPS 

SDsR’C^t^SO 

Ci'1  = S0*CM 

SN=-OM*EXP2 

GO  TO  4l 

42  K=K/2 

43  continue 

C OIVIDES  EACH  ELEMENT  OF  X2  8Y  N WHEN  TRANSFORM  IS  FROM  TIME 
C SEQUENCE  TO  FREQUENCY  SPECTRUM. 

50  ifcisgno.gt.o)  return 

CF=1.0/FLOAT(Ni ) 

00  51  K=1,NT2 

X2(K) =CF*X2(K) 

51  CONTINUE 
RETURN 

C ERROR  EXIT 

90  WRITE  (LUN,2090)  NO 
RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 
2090  FORMAT I1X/24H  IMPROPER  N VALUE. /faH  N =,15/ 

1 35H  ERROR  DETECTED  IN  ROUTINE  OFT  /) 

END 


SUBROUTINE  FRBPFC (NOB,RP8,FC0,eWB.UFS. IFCMNtlFCMX.HFB) 

C THIS  SUBROUTINE  CALCULATES  FREQUENCY  RESPONSE  OF  A BANDPASS 
C FILTER  (BPF)  OF  CHEBYSHEV  TYPE.  (CF.  STORER.  PASSIVE  NETWORK 
C synthesis,  MCGRAW-HILL,  1957) 

c the  input  parameters  are 

C NOB  = ORDER  OF  THE  BPF 

C (MUST  NOT  BE  NEGATIVE  NOR  GREATER  THAN  100), 

C RPB  = RIPPLE  (03)  WITHIN  THE  PASSBAND  OF  THE  BPF 

C (MUST  BE  NONNEGATIVE) , 

C FCB  = CENTER  FREQUENCY  OF  THE  BPF 

C (MOST  BE  POSITIVE!  , 

C BW8  = 3-DB  BANOWIOTH  OF  THE  BPF  (MUST  BE  POSITIVE), 

C UFS  s UNIT  FREQUENCY  SPACING  (MUST  BE  POSITIVE), 

C IFCMN  = MIN.  OF  FREQ.  COMPONENT  NUMBER, 

IFCMX  * MAX.  OF  FREQ.  COMPONENT  NUMBER. 


C THE  OUTPUT  parameter  IS 

C HF3  = complex  array  FOR  THE  FREQUENCY  RESPONSE  OF 
C THE  BPF. 

C DECLARATION  STATEMENTS 

complex  HFB(IOO) 

complex  pole (100) ,HF3I,HFC,NMRTR,Z 
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DATA  NOdMX/lOC/ 

DATA  NOPV/-1/,  RPPV/-1.0/ 

DATA  HLFPI/1.3/j79  63267/ 

DATA  DRMN/l.iic-ia/ 

DATA  LUN/b/ 

ASINH<X>=ALO6(X»SURT{X*X*1.0)  » 

C PRELIMINARY  PROCESSING 
■ IG  N0=N03  , 

RP=RP8 

FC=FC3 

3M=dWd 

UF=UFS 

IMN=IFCMN 

IMX=IFCMX 

IFINO.LT.O.OR.NO.Gl.NOeMXJ  GO  TO  90 
IF(RP.LT.O.L)  go  TO  9ii 

IF (FC.LE.O.C)  GO  TO  90 

IFOW.LE.O.O)  GO  TO  90 

IFIUF.LE.O.O)  GO  TO  90 

IF(NO.EU.O)  GO  TO  70 

C CALCULATES  BASIC  PAkAMETERS  HHEN  NECESSARY. 

20  IFINO.EQ.NOPV. ANO.RP.EO.RPPtf)  GO  TO  60 

C - NOB.GE.3. AND.RPQ.GT.0.0 

NOPO=NO 

RPPrf=RP 

NOPl=NO»l 

FLNO=NO 

R=1.0 

HFC=1.0 

PHI0=HLFPI/FLNO 

IFCRP.EQ.O.O.OR.NO.EQ.lt  GO  TO  40 

VR=EXP<0.11512925>.6*RP) 

RE  = 1.0/SQRT  TV«*y/K-1.3) 

U0=ASINHTRE)/FLNO 
EXU=EXP(U0t 
REXU=1. 0/EXU 
CH=0.5*TEXU*REXUI 
SH=0.5* (EXU-REXU) 

IF (N0-N0/2*2.EG,0)  HFC=HFC/VR 
lF(NO.EQ.2t  GO  TO  30 

21  TH0R=R»R 
TM2  = R 

TM1=THOR*R-1.0 


OM2=1.0 
OM1=TMOR»THOk 
DO  22  K=3.N0 

T=THOR*TMl-TM2 
0aTMl»TMl»TW0R*0Ml-0M2 
TM2=TM1 
TM1=T 
OM2=OMl 
0Ml=O 

22  CONTINUE 
RPV=R 

R=R-  (T-RE)  /£) 

IF (ABS(R-RPV) .GT.ORMN)  GO  TO  21 
GO  TO  50 

C - NOB.EQ.2 

30  R:SQRT(0.5*(1.04RE> ) 

GO  TO  50 

C * NOB . EQ • 1 . OR . RPB • E Q • 0 « 0 
40  CM=1.0 
SHsl.O 

C CALCULATES  POLES. 
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I 

I 


— * - 


50  nmrtr=hfc 

00  5X  K=1,N0 

PHI=FL0AT (KtK-NJPl) *PHI 0 

POLE IK) =CMPLXI -SH*C0S IP  HI) ,CH*5INIPm1)  ) 

NMRTR=-NMRTR*POLE IK) 

51  CONTINUE 

c calculates  frequency  response. 

. 60  CF=2. a*R/aW 

do  62  I=IMN)IHX 

Z=CHPLXI O.U , IFLOAT I I-l) »UF-FC) »CF) 

HFBI=NMRTR 
00  61  K=1,NC 

MF8I  = HFBI/ (2-POLE  IK) ) 

61  CONTINUE 
HFOIDsHFBI 

62  CONTINUE 
RETURN 

c «LL  PAj;  filter  INOa=0> 

•70  00  71  I = IMN,IMX 
HF8(I)=I1. 0,0.0) 

71  CONTINUE 
RETURN 
C ERROR  EXIT 

90  WRITE  (LUN,209Q)  NO , RP, FC. 3M , UF 
RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 
2090  FORMAT  I1X/31H  »*♦  IMPROPER  INPUT  ylALUEIS)./ 

1 3X,5HN03  = , I5,5X,5HRPB  = , £ 11 . 3 , 5 X, 5HF CB  =,E11.3, 

2 5X,5H3WB  =,E11. 3,5X,5HUFS  =,E11.3/ 

3 35H  ERROR  OtTECTEO  IN  ROUTICE  FRBPFC/) 

END 


SUBROUTINE  MOOAM  I NF C , UFS , SF  ,FC , FLMI , VF ) 

THIS  SUBROUTINE  SIMULATES  AN  AMPLITUDE-MODULATION  lAM) 
LATOR. 

THE  INPUT  PARAMETERS  ARE 


MOOU- 


NFC 

UFS 

SF 


FC 

FLMI 


NUMBER  OF  FREDHEnCY  COMPONENTS, 

UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 

COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 
tfOLTAGE  SPECTRUM  COMPONENTS  OF  THE  MODULATING 
SIGNAL, 

CENTER  FREQUENCY  OF  THE  MODULATED  SIGNAL, 
MODULATION  INDEX. 

THE  OUTPUT  PARAMETER  IS 

VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  MODULATED  SIGNAL  ARE 
TO  BE  STORED. 

THIS  subroutine  CALLS  THE  RANF  FUNCTION. 

DECLARATION  STATEMENTS 

DIMENSION  SF  (6<»)  , VFI64) 

REAL  MI 

DATA  SQRT2/1.<,1<,2135623  7310/,TWOPI/6.28  318  530  71795  9/ 
DATA  LUN/6/ 

PRELIMINARY  PROCESSING 
10  N0=NFC 
U0=UFS 


F0  = FC 
MI=FLMI 


IF  IN0.lt. <,) 

GO 

TO 

90 

IFIUO.LE.O.C) 

GO 

TO 

9C 

IFC=F0/Ui*1.5 

IFIIFC.GE.NO) 

GO 

TO 

90 

1 


i 

I 

! 


1 

j 
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IF(MI.&T.1.C)  00  TO  40 

NT2=NQ*2 

IFCT2=IFC*2 


C GFNtRATES  THE  CARRIER  PHASE. 

20  PHI  = THOPI*<R.ANF(0)-0.5) 

CP=CCS«PHI) 

sp=s:n(phu 

CCETERMINES  The  COEFFICIENTS  FOR-  SIOEBANO  SPECTRAL  COMPONENTS. 
30  PS=a.3 

00  31  J=3,NT2 

PS=PS*SF(J)**2 
. 31  CONTINUE 

PS=PS/2. 0 

IF(PS.LE.G.O)  GO  TO  91 
C=0.5*MI/SCRT(PS> 

C1=C*CP 

C2=C*SP 

C calculates  EACH  SPECTRAL  COMPONENT  OF  THE  MODULATED  SIGNAL 
C and  STORES  THE  RESULTS  IN  THE  \IF  ARRAY. 

C - DC  COMPONENT. 

AO  VF(1)=0.0 
y/F(2)  = 3.; 

C - LOWER  SIDEBAND. 

A2  KMX=IFCT2-2 
J=IFCT2 

DO  43  K=4.KMX|2 

J = J-2 

S1=SF<J-1) 

S2  = SF  « J» 

>/FIK-l)=  C1*S1+C2»S2 
VF(K1  =-Cl*S2+C2*Sl 

43  CONTINUE 
C - CARRIER. 

44  VF (IFCT2-1)=CP*SURT2 
\/F(IFCT2»  =SP*SQRT2 

C - UPPER  SIDEBAND. 

46  KMN=IFCT2f2 
J = 2 

00  47  K=KMN.NT2,2 
J=J*2 
S1=SF ( J-1) 

S2=SFC J) 

VF(K-1)=  C1*S1-C2*S2 
VF(K)  = C1*S2*C2*S1 

47  CONTINUE 

return 

C ERROR  EXIT 

90  WRITE  (LUN,2090)  NO.UO.FO 
GO  TO  92 

91  WRITE  (LUN,Z091) 

92  WRITE  CLUN,2092) 

return 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 

2090  F0RMAT<1X/31H  IMPROPER  INPUT  VALUE(S»./ 

1 8H  NFC  =.I6.5X,5HUFS  =. Ell . 3 , 5X, 4HFC  =,E11.3) 

2091  F0RMAT(1X/33H  ***  ALL  ZERO  COMPONENTS  IN  SF.) 

2092  FORMAT(35H  ERROR  DETECTED  IN  ROUTINE  MOOAM  /» 

ENO 


SUBROUTINE  MOOFM  { NFC , UFS , SF  ,FC ,FO V , VF) 

C This  subroutine  simulates  a frequency-modulation  (FM)  moou- 

C LATOR. 
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C THE  INPUT  PARAMETEfS  ARE 


NFC  = NUMBER  OF  FREQUENCY  COMPJNENTS, 

UFS  = UNIT  FRECUENCY  SPACING  OF  THE  SPECTRUM, 

SF  = COMPLEX  ARRAY  OF  OIMENSION  NFC  CONTAINING  THE 
VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  MOOULATING 
SIGNAL, 

FC  = CENTER  FREQUENCY  OF  THE  MODULATED  SIGNAL, 

FDV  = FREQUENCY  DEVIATION.  . 

THE  OUTPUT  PARAMETER  IS 

VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  MODULATED  SIGNAL  ARE 
TO  BE  STORED. 

THIS  SUBROUTINE  CALLS  THE  OFT  SUBROUTINE  AND  THE  RANF 
FUNCTION. 

DECLARATION  STATEMENTS 

DIMENSION  SF(6<»)  ,VF(q<» 

DATA  SURT 2/ 1.4 11.2 13 562 3 73 lu/,TrtOPI/6. 28  3185  30  71795  9/ 
DATA  LUN/6/ 

C PRELIMINARY  PROCESSING 
10  N0=NFC 
U0=UFS 
FO-FC 
FO=FOV 


IF(N0.LT.4) 

GO 

TO 

90 

IF(UO.LE.O.C) 

GO 

TO 

90 

IFC=F0/U0»1.5 

IFdFC.GE.NC) 

GO 

TO 

90 

FO  = UO*FLO.*T(IFC-1) 

NT2=N0*2 

DOMG=TWOPI*Uu 

C INTEGRATES  THE  MOOULATING  SIGNAL  IN  THE  FREQUENCY  DOMAIN, 

C STORES  THE  RESULTS  IN  THE  VF  ARRAY,  AND  DETERMINES  THE 
C COEFFICIENTS  FOR  PHASE  MODULATION. 

30  VF(1)=0,J 
VF«2)=0.0 
PS  = 0. 0 
OMG=0.0 

DO  31  J=4,NT2,2 

OMG=OMG+OOMG 
VF(J-1)=  SF(J)  /OMG 
VF(J)  =-SF(  J-D/OMG 
PS=PS»SF  <J-1)**2*SF (J)**2 

31  CONTINUE 

IF (PS.LE.O.OI  GO  TO  91 
C=FO/SQRTIPS) 

C FOURIER-TRANSFORMS  THE  VF  ARRAY  FROM  FREQUENCY  DOMAIN  TO  TIME 
C DOMAIN. 

40  CALL  OFT(NO,VF,1,VF> 

C PHASE-MODULATES  THE  CARRIER,  AT  EACH  SAMPLING  POINT  IN  THE 
C TIME  DOMAIN,  WITH  THE  VF  ARRAY  AND  STORES  THE  RESULTS  IN  THE 
C VF  ARRAY. 

50  FMX=U0*FLOAT<N0) 

CC=RANF«0) 

FI=-F0 

DO  51  J»2,NT2,2 

FI=FI*F0 
CYCL*FI/FMX*CC 
CYCL=CYCL-AINT (CYCLIC. 51 
PHI  = TWOPI*«CYCL«-C*VF(  J-ll  ) 

VF( J-l)=SQRT2*COS(PHI) 

VF<J)  =SQKT2*SINIPHU 

51  CONTINUE 

C FOURIER-TRANSFORMS  THE  MODULATED  SIGNAL  FROM  TIME  DOMAIN  TO 
C FREQUENCY  DOMAIN. 
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60  CALL  OFT(NOfVFt-lt\/F» 

RETURN 
C ERROR  EXIT 

90  WRITE  (LUN,E09o)  NJtUOtFO 
GO  TO  92 

91  WRITE  (LUN,2091) 

92  WRITE  ILUN,2092) 

RETURN 

C FORMAT  STATEMENT'S  FOR  ERROR  MESSAGES 

2090  FORMAT ( 1X/31H  ***  IMPROPER  INPUT  VALUE  (S)./ 

1 6H  NFC  =,Ib,aX,5HUFS  = , £11 . 3 , 5Xt  <*HFC  =fE11.3) 

2091  FORMAT (1X/33H  ♦**  ALL  ZERO  COMPONENTS  IN  SF.) 

2092  FORMAT(35H  ERROR  DETECTED  IN  ROUTINE  MOOFM  /) 

END 


SUBROUTINE  MODFSK(NFC,UFSfNlB,IBfFC.FS, VF) 

THIS  SUBROUTINE  SIMULATES  A FRLQUENCY-SHIFT-KEYING  (FSKI  MODU- 
LATOR. 

THE  INPUT  PARAMETERS  ARE 

NFC  = NUMBER  OF  FREOUENCY  COMPONENTS, 

UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 

NIB  = NUMBER  OF  INFORMATION  BITS  IN  A PERIOD 
(MUST  NOT  EXCEED  NFC/21, 

IB  = ARRAY  OF  DIMENSION  NIB  CONTAINING  THE  INFORMA- 
TION BITS, 

FC  = CENTER  FREQUENCY  OF  THE  MODULATED  SIGNAL, 

FS  = FREQUENCY  SHIFT. 

THE  OUTPUT  PARAMETER  IS 

VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  MODULATED  SIGNAL  ARE 
TO  BE  STORED. 

THIS  SUBROUTINE  CALLS  THE  OFT  SUBROUTINE  AND  THE  RANF 
FUNCTION. 

DECLARATION  STATEMENTS 

DIMENSION  13 (16) , VF 164) 

DATA  SQRT2/1.41*,21356237310/,  TWOPI/6. 28318530717959/ 

DATA  LUN/b/ 

C PRELIMINARY  PROCESSING 
13  N0=NFC 
U0=UFS 
NI=NIB 
FCO=FC 
FS0=FS 

IF(NU.LT.4)  GO  TO  90 

IF(UC.LE.O.O)  GO  TO  90 

IF(2*NI.GT.N0)  GO  TO  90 

FO=UO*4INT(FCO/UO+3.5)-0.5*FSO 

F1=F0*FS0 

FMX=U0*FLOAT(NC) 

IF(FO.LE.O.O)  GO  TO  90 
IF (Fl.GE.FMX)  GO  TO  90 
FLNI=NI 
OLT=1.0/FHX 

C FREQUENCY-SHIFT-KEYS  THE  CARRIER,  AT  EACH  SAMPLING  POINT  IN 
C THE  TIME  DOMAIN,  BY  CHANGING  THE  CARRIER  PHASE  ACCORDING  TO 
C THE  INFORMATION  BIT  SEQUENCE  AND  STORES  THE  RESULTS  IN  THE 
C VF  ARRAY. 

50  CYCL=RANF(0» 

I8IPV=IB(1) 

DO  59  1=1, NO 

IF(I.EQ.l)  GO  TO  58 

IIB=( ( I-1»*NI)/N0+1 


orjoooooooooooooo 


f 


iai=i3(iiB) 

IFtI9I.£U.I3IPV)  GO  TO  54 
51  I3IPy/=IBI 

TT=FLO«T  (NO*  ( II3-in  / FLNI 
0LT1=  tTT-FLOAT<I  -2n*0LT 
0LT2=0LT-DLT1 

IFdai.GT.C)  GO  TO  53 

• 52  CYCL=CTCL*F1*OLT1»FO*OLT2, 

GO  TO  5iJ 

53  CYCL=CYCi.*Fa*OLTl*Fl*OLT2 
GO  TO  58 

54  IF(  IBI.GT.O)  GO  TC  5t 

55  CYCL=CYCL*FO*OLT 
GO  TO  58 

56  CYCL=CYCL*F1*DLT 

58  CYCL=CYCL-AINT (CYCL *0 . 51 
PHI=THOPI*CYCL 

J = I*I 

YF  IJ-l)=SQRT2*COSIPHn 
VF(J)  sSQRT2*SINlPHI) 

59  CONTINUE 

C F0URIER-TRANSF0RM3  THE  MOOULATtO  SIGNAL  FROM  TIME  DOMAIN  TO 
C FREQUENCY  DOMAIN. 

60  CALL  OFT(NO.VF,-1, YF) 

RETURN 

C ERROR  EXIT 

90  MRITE  (LUN, 20901  NO . U J , NI , FC o «F SO 
92  WRITE  (LUN.2092) 

RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 
2090  FORMAT (1X/31H  ***  IMPROPER  INPUT  YALUE(S)./ 

1 8H  NFC  =, Ifa,5X,5HUFS  = , E 1 1 . 3 .5 X , 5HN 1 B =,Ib, 

2 5X.5HFC  =»cll.3»5X,5HFS  =,£11.3) 

2092  FORMAT(35H  ERROR  DETECTED  IN  ROUTINE  MOOFSK/) 

END 


SUBROUTINE  M0LPHM(NFC,UFS,SF,FC,FLMI,YF) 

THIS  SUBROUTINE  SIMULATES  A PHASE-MODULATION  MODULATOR. 

The  INPUT  PARAMETERS  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS, 

UFS  = UNIT  frequency  SPACING  OF  THE  SPECTRUM, 

SF  = complex  array  of  DIMENSION  NFC  CONTAINING  THE 
VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  MODULATING 
SIGNAL, 

FC  = CENTER  FREQUENCY  OF  THE  MODULATED  SIGNAL, 

FLMI  = MODULATION  INDEX  (IN  RADIANS). 

THE  OUTPUT  PARAMETER  IS 

YF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  MODULATED  SIGNAL  ARE 
TO  BE  STORED. 

THIS  SUBROUTINE  CALLS  THE  OFT  SUBROUTINE  AND  THE  RANF 
FUNCTION. 

DECLARATION  STATEMENTS 

DIMENSION  SF(b4)  ,YF(t>4) 

REAL  MI 

DATA  SQRT2/1.4142135623731u/,  T WOPI /6. 2 8 318530 7 17959/ 
DATA  LUN/6/ 

C PRELIMINARY  PROCESSING 
10  N0=NFC 
UC=UFS 
F0  = FC 
MI=FLMI 


J 


1 

I 

! 

i 

\ 


i 

I 

I 

i 


I 

i 

i 

\ 

i 

! 

1 


50 


I 

] 
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IFIN0.lt. h) 

GO 

ro 

90 

IF(UC.L£.U.0> 

GO 

TO 

90 

IFC  = F0/U0H.5 

IFtlFC.GE.NO) 

GO 

TO 

93 

FO=UO*FLOAT (IFC-l) 


NT2=N0*2 

C OETEkMINES  THE  COEFFICIENT  FOR  PHASE  MOOULATION. 

30PS  = ij>J 

00  31  J=1,KT2 

PS=PS*SF(J)**2 
31  CONTINUE 

IFJPS.LE.t.O)  GO  TO  91 
C=MI/SaRT(PS) 

C FOURIER-TRANSFORMS  THE  MCOUCATINU  SIGNAL  FROH  FREQUENCY  DOMAIN 
C TO  TIME  DOMAIN  AND  STORES  THE  RESULT  IN  THE  VF  ARRAY. 

43  CALL  OFTCNO.SF.IWF) 

C PHASE-MODULATES  THE  CARRIER,  AT  EACH  SAMPLING  POINT  IN  THE 
C time  DOMAIN,  WITH  THE  VF  ARRAY  ANU  STORES  THE  RESULTS  IN 
C the  VF  ARRAY. 

50  FMX=UJ*FLOAT(N3) 

CC=RANF(-1) 

FI=-F0 


00  51  J=2,NT2,2 

FI=FI»FO 
CYCL=FI/FMX+CC 
CYCL=CYCL-AINT (CYCLED. 5) 

PHI=TW0PI*CYCL*C*VF(J-1) 

VF< J-1)=SURT2*C0SCPHI) 

VF(J)  =SQRT2*SIN(PHI) 

51  CONTINUE 

C FOURIER-TRANSFORMS  THE  MODULATED  SIGNAL  FROM  TIME  DOMAIN  TO 
C FREQUENCY  DOMAIN. 

60  CALL  OFT(Na,VF,-l,VF» 

RETURN 
C ERROR  EXIT 

90  WRITE  (LUN,209u)  NO,UO,FO 
GO  TO  92 

91  WRITE  (LUN,2091» 

92  WRITE  (LUN,2092» 

RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 

2090  FORMATtlX/31M  **♦  IMPROPER  INPUT  VALUE(S)./ 

1 3H  NFC  =,I6,5X,5HUFS  = , Ell . 3 , 5X , 4HFC  s,Ell.3> 

2091  F0RMAT(1X/33H  *•*  ALL  2ERO  COMPONENTS  IN  SF.) 

2392  FORMAT(35H  ERROR  DETECTED  IN  ROUTINE  MOOPHM/J 

END 


SUBROUT  INF.  MOOSSiHNFC,UFS,SF,FR,  VF) 

THIS  SUBROUTINE  SIMULATES  A SINGLE-SIDEBAND  (SSB)  MODULATOR, 
THE  INPUT  PARAMETERS  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS, 

UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 

SF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 
VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  MODULATING 
SIGNAL, 

FR  = REFERENCE  FREQUENCY  OF  THE  MOOJIATEO  SIGNAL. 

THE  OUTPUT  PARAMETER  IS 

VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  MODULATED  SIGNAL  ARE 
TO  BE  STORED. 

THIS  SUBROUTINE  CALLS  THE  RANF  FUNCTION. 

DECLARATION  STATEMENTS 
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DIMENSION  SF(6<*)  , \/F(64) 

DATA  TW0PI/6.i.’93H53  J717'J59/ 
DATA  LUN/b/ 


C PRELIMINARY  PROCESSING 
10  N0=NFC 
UC=UFS 
F0  = FR 

IF(NC.LT.<.)  GO  TO  90  , 

IF(UO.LE.O.O)  GO  TO  90 

IFR=FQ/U0+1.5 

IF(IFR.GE.NC)  GO  TO  90 

NTE=N0*2 

IFRT2=IFR*2 

C GENERATES  THE  CARRIER  PHASE. 

20  PHI=TWOPI*<RANF 10)-3.6> 

CP=C05TPHI) 

SP=SIN(PHI» 

C OETERMINES  Trit  COEFFICIENTS  FOR  iSb  HOOOLATICN. 

•30  PS  = 0.0 

DO  31  J=3,NT2 

PS=PS»SF«J)**2 
31  CONTINUE 
PS=PS/2.0 

IF  (PS. L£. 0.0)  GO  TO  91 
C=1.0/SQRT(FS) 

C1=C*CP 

C2=C*SP 

C TRANSLATES  EACH  COMPONENT  OF  THE  MODULATING  6ASESAN0  SIGNAL 
C UP  TO  THE  RF  BAND  AND  STORES  THE  RESULTS  IN  THE  y/F  ARRAY. 

<♦0  DO  41  K = 1,1FRT2 
\/F(K)  = 0.0 

41  CONTINUE 

46  KMN=IFRT2»2 
J = 2 

DO  47  K=KMN,NT2,2 
J = J + 2 

S1=SF(J-1) 

S2=SF( J) 

\/F<K-l»=Cl*Sl-C2*S2 
VF(K)  =C1*S2*C2*S1 

47  CONTINUE 
RETURN 

C ERROR  EXIT 

90  WRITE  (LUN,2D9U  NO.UO.FO 
GO  TO  92 

91  WRITE  (LUN,2C91) 

92  WRITE  tLUN.2092) 

RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 

2090  F0RMAT(1X/31H  •**  IMPROPER  INPUT  \/ALUE(S>./ 

1 8H  NFC  =.Ifa.5X,5HUFS  = , E 11 . 3, 5X, 4HFR  =,E11.3) 

2091  FORMAT (1X/33H  *•*  ALL  ZERO  COMPONENTS  IN  SF.) 

2092  F0RMATC35H  ERROR  OETtCTEC  IN  ROUTINE  MODSSB/) 

END 


SUBROUTINE  PSPHCT < NFC . UFS . SF ,PW ) 

C this  subroutine  SIMULATES  A PSOPHOMETER  FOR  COMMERCIAL- 
C TELEPHONE  CIRCUIT.  (CF.  CCITT  RECOMMENDATION  P.53) 

C THE  INPUT  PARAMETERS  ARE 

C NFC  = NUMBER  OF  FREQUENCY  COMPONENTS, 

C UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 

C SF  = COMPLEX  ARRAY  OF  QIMENSIUN  NFC  CONTAINING  THE 


C y/OLTA&E  SPECTRUM  COMPONENTS  OF  THE  INPUT  SIGNAL. 

C ThE  OUTPUT  PARAMETER  IS 

C PW  = PSOPHOMETRICALLY-HEIGHTEO  SIGNAL  POWER. 

C DECLARATION  STATEMENTS 


DIMENSION 

SF(64) 

OIMENSION 

WT<500I 

DIMENSION 

ALFI30I , 

08(301 

DATA  ALF 

.■ 

1 /I. 22185. 

1.69897, 

2. 00000, 

2.17609, 

2.30103, 

2 2.39794, 

2.47712, 

2.54407, 

2.60206, 

2.65321. 

3.  2.69897, 

2.77815, 

2.84510. 

2.90309, 

2.95424, 

4 3.00000. 

3.04139, 

3.07913, 

3.11394, 

3.17609, 

5 3.30103, 

3. 34242, 

3.39794, 

3.47712, 

3.54407, 

6 3.56620, 

3.60206, 

3.63347, 

3.65321, 

3.69897/ 

DATA  08 

1 / >85.00, 

-63.00, 

-41.00, 

-29.00 , 

-21.00, 

2 -15.00, 

-10.60, 

-8.50. 

—6.30, 

-4.70, 

3 —3.60, 

-2.00, 

-0.90, 

0.00, 

0.60, 

4 1.00, 

0.60, 

0.00, 

-0.40, 

-1.30, 

5 -3.00, 

- 3 . 8 f 

-4.20, 

-5.60, 

-8.50, 

6 -10.70, 

-15.0  0, 

-20.70, 

-25.00, 

-36.30/ 

DATA  NPV/O 

/,  UPy/70. 

0/,  FMAX/IOOOO.O/, 

NMAX/500/ 

C PRELIMINARY  PROCESSING 

10  N0=NFC-1 
U0=UFS 

FMX=UO*FLOATINC I 
IF(FMX.LE.FMAX)  GO  TO  11 

FMX=FMAX 
N0=FMX/U0 

11  NG=1 

15  IF(ND.LE.NMAX)  GO  TO  20 

N0=(N0-l>/3 
NG=NG*3 
UO=UQ»3.0 
GC  TO  15 

C CALCULATES  THE  WEIGHTING  COEFFICIENTS  WHEN  NECESSARY. 
20  IFINO.LE.NPV.ANO.UO.ECl.UPV)  GO  TO  30 
NPW=N0 
UPV=U0 
JPV  = 0 


J=2 

FRI=0.0 
00  29  I=1.N0 

FRI=FRI+U0 

IFIFRI.LE.C.COOl)  GO  TO  21 

ALFI=0.43*»29<.i,819*ALOG»FRI) 

GO  TO  22 

21  ALFI  = -«».0 

22  IFIALFI.lt. ALFIJ) ) GO  TO  23 

IF(J.GE.3b)  GO  TO  23 

J=J*1 
GO  TO  22 

23  IFIJ.EQ.JPV)  GO  TO  24 

JPtf=J 

A0=ALF<J-1) 

C0s08«J-l) 

Cl::  l09IJ)-0dlJ-l)  l/(ALFf  Ji-ALFI  J-ll  ) { 

24  WTI=C0»C1*<ALFI-A0)  j 

WT(I)=EXP<0.23a2585093«UTIl  f 

29  CONTINUE 

C CALCULATES  THE  WEIGHTEO  POWER.  I 

30  P0=0.0 
IFCsNG 

DO  39  I<ltN0 
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j IFC=IFC»a 

P0=P0»WT<II*(SFaFC»**2»SF(IFC*ll**2l 
^ 3d  CONTINUE 

39  CONTINUE 

PW=J.5*P0*UFS 

: RETURN 

t . END 


SUBROUTINE  PSPHPT<NFC,UFS,SF,PH) 

THIS  SUBROUTINE  SIMULATES  A PSOPHOMETER  FOR  PROGRAM- 
TRANSMISSION  CIRCUIT.  TCF.  CCITT  RECOMMENDATION  P.53) 

THE  INPUT  parameters  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS* 

UFS  = UNIT  FREQUENCY  SPACING  OF  ThE  SPECTRUM, 

SF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  CONTAINING  THE 

YOLTAGE  SPECTRUM  COMPONENTS  OF  THE  INPUT  SIGNAL. 
The  OUTPUT  PARAMETER  IS 

pw  = psophometrically-w.'ighted  signal  power. 

DECLARATION  STATEMENTS 
dimension  SFJ6L) 

DIMENSION  HT(500)- 
dimension  ALF(30) ,09(30) • 

DATA  ALF 


1 

/l.GJOOJ, 

1.30103, 

l.SgTBi,, 

1.47712, 

1.60206, 

z 

1.69897, 

1.77815, 

1.90309, 

2.00000, 

2.17609, 

3 

2.30103. 

2. *,7712, 

2.60206, 

2.77813 , 

2.90309, 

h 

3.00000. 

3.17609, 

3.30103, 

3.47712, 

3.60206, 

5 

3.69897, 

3.77815, 

3. 84510 , 

3.90309, 

3.95424, 

6 

<,.00000, 

0.11391,, 

4.176C9, 

4.30103, 

4.60206/ 

DATA  DB 

1 

/ ->,1.50, 

-1.1.50, 

-41.00, 

-40.00, 

-37.50, 

Z 

-3>,.30, 

-32.20, 

-29.10, 

-26.10, 

-20.60, 

3 

-17.30, 

-12.70, 

~ d • 6 0 9 

-4.80, 

-1.90, 

L 

0.30, 

3.20, 

5.30t 

7.00  , 

8.20, 

5 

3.<,0  , 

8.20, 

7.30« 

5.10, 

-0,30, 

6 

-9.70, 

* 3 0 • Q D t 

* 3 3 « 0 li  f 

-36.00, 

-40.00/ 

DATA  NPY/a/ 

, UPJ/O.i 

0/,  FMAX/hOCuO.O/, 

NMAX/500/ 

C PRELIMINARY  PROCESSING 

10  N0=NFC-1 
U0=UFS 

FMX=UJ*FLOAT (NO ) 

IF(FMX.LE.FMAX)  go  to  11 

FMX=FMAX 

N0=FMX/U0 

11  NG=1 

15  ifjnc.le.nmaxi  go  to  ZO 

N0=(N0-l)/3 
NG=NG*3 
U0=U0*3. 0 
GO  TO  15 

C CALCULATES  THE  WEIGHTING  COEFFICIENTS  WHEN  NECESSARY, 
ZO  IF(N0.LE.NPV.AND.U0.£Q.UPY)  go  TO  30 
NPV=N0 
UPY=U0 
JPV  = 0 
J=Z 

FRIsO.O 
00  Z9  1=1, NO 
FRI=FRI*U0 

IF (FRI.LE.C. 0001 ) GO  TO  Z1 

ALFI=0,4  3i»2  944  819*ALOG(FRI) 
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liO  TO  22 

21  ALFI=-i,.iJ 

22  IF(AlFI.LT. ALFIJ) ) GO  TO  23 

IF(J.GE.3&)  GO  TO  23 

J = J*1 
GO  TO  22  • 

23  IFtJ.EU.JPV)  GO  TO  2A 

JP\I  = J 

A0=ALF(J-1) 

C&=03(J-1) 

Cl= 1091 J)-D3 (J-l) )/ (ALF IJ)-ALF  « J-l»  » 

Z>*  WTI=C3*C1*<ALFI-A0) 

WT<II=EXP(0.23C2585093*MTI) 

29  CONTINUE 

C CAUCULATtS  THE  WEIGHTED  POWER. 

30  PG=0.3 
IFC=NG 

DO  39  1=1. NO 

DO  39  J=1.NG 
IFC=IFC*2 

P0  = P0»WT<U*<5F<IFCJ**2»SF(IFC*1I**2» 
39  continue 
39  CONTINUE 

PW=0.5*PJ*UFS 

RETURN 

END 


SUBROUTINE  SGPLS  < NFC, UFS ,FC ,FRF ,TR, TW , TF , 0 VF, VF) 

THIS  SUBROUTINE  SIMULATES  A PULSE  jIGNAL  GENERATOR. 

THE  INPUT  PARAMETERS  ARE 

NFC  = NUMBER  OF  FREQUENCY  COMPONENTS. 

UFS  = UNIT  FREQUENCY  SPACING  OF  THE  SPECTRUM, 

FC  = CENTER  FREQUENCY  OF  OUTPUT  SIGNAL, 

PRF  = PULSE  REPETITION  FREQUENCY, 

TR  = PULSE  RISE  TIME, 

TW  = PULSE  WIDTH  AT  THE  TOP, 

TF  = PULSE  FALL  TIME, 

DVF  = MAXIMUM  FREQUENCY  OEYIATION, 

THE  OUTPUT  PARAMETER  IS 

\IF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE  VOLTAGE 
SPECTRUM  COMPONENTS  OF  THE  OUTPUT  SIGNAL  ARE  TO 
3E  STORED. 

THIS  SUBRDUTINE  CALLS  THE  DFT  SUBROUTINE  AND  THE  RANF 
FUNCTIDN. 

DECLARATION  STATEMENTS 
DIMENSION  VF<G4) 

DATA  SQRT2/1.41421356237310/,TWOPI/6. 28 31 853071795 9/ 
DATA  LUN/6/ 

C PRELIMINARY  PROCESSING 
10  N0=NFC 
U0»UFS 
F0  = FC 
P0=PRF 

IF(N0.LT.4)  GO  TO  90 

IFiUO.LE.ti.Oi  GO  TO  90 

IF  IPG. LE. 0.0)  GO  TO  90 

TROsTR 

TW0=TW 

TF0=TF 

00=OVF 

FMX=UO*FLOAT(NO» 

P0=U0«AINT (P0/U0»0.5) 


1 

( 

t ‘ 

I 


F8=F0-JJ 

Ft=F0»3O 

IF»F8.U£.0.0.0(<.F3.G£.FMXI  GO  TO  90 

if(Fl.le.o.o.ok.fe.ge.fnx*  go  to  90 

PRO=1.3/P0 

T2=JPRO-TwO)/2.C 

T1=T2-TR0 

T3=T2+Trt0 

T4=T3fTF0 

IF(Tl.LT.O.O)  GO  TO  9u 

IF  (T<..GT.PRC)  GO  TO  90 

T4i  = T-*-Tl 

IF(T1*1.LE.0.0)  GO  TO  9u 

C DETERMINES  THE  COEFFICIENTS  FOR  PHASE  MODULATION. 

2U  C0=RANF(0) 

C1=(F3*TA-FE*T1)/T<.1 

C2=(FE-F0)/«2.C*T41) 

C GENERATES  THE  PULSE  SIGNAL  AT  EACH  SAMPLING  POINT  IN  THE  TIME 
C DOMAIN  ANO  STORES  THE  RESULTS  IN  THE  t/F  ARRAY, 

30  CO  39  1=1, NO 

J = I+I 

T = FLOATlI-n/FiiX 
Tfl=AMOOIT ,PROJ 
IF(TO.LE.TI)  go  TO  38 

IF(T0.GE.T4»  GO  TO  '36 

IF  (TJ.GT.T3)  GO  TO  33 

IF(T3.GE.T2)  GO  TO  32 

31  A=(TC-T1)/TR0 
GO  TO  34 

32  A=1.0 

GO  TO  34 

33  A=U>*-T0>/TFC 

34  A=SQRT2*A 
CYCL=C0FT*IC1+T»C2) 

CYCL=AMOOICYCL+0.5, 1.01-0.5 
PHI=TWOPI*CYCL 
VF<J-l)=A*COS (PMl> 

VF(J)  =A*SIN(PHI) 

GO  TO  39 

38  VF(J-1)=0.C 
OF(J)  =0.0 

39  CONTINUE 

C FOURIER-TRANSFORMS  the  generated  signal  from  TIME  DOMAIN  TO 
C FREQUENCY  DOMAIN. 

40  CALL  OFT(NO,VF,-1, 0F> 

RETURN 

C ERROR  EXIT 

90  HRITE  CLUN.2090)  NO , UO , FO , PC , TRO , TNO , TF 0 • DO 
RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 
2090  F0RMAT(1X/31M  IMPROPER  INPUT  VALUE <S)./ 

1 3X.5HNFC  =,Ifc,  5X.5HUFS  =, £ 1 1. 3, 5X, 5HFC  =,E11.3, 

2 5X,5HPRF  =.E11. 3/3X,5HTR  =, El 1. 3, 5X , 5HT H =,E11.3, 

3 5X.5HTF  =,E11. 3,5X,5H0VF  =,E11.3/ 

4 35H  ERROR  DETECTED  IN  ROUTINE  SGPLS  /) 

END 


SUBROUTINE  SGWGN  ( NFC, UFS , IF CMN, IFCMX. VF 1 
C THIS  SUBROUTINE  SIMULATES  A WHITE  GAUSSIAN  NOISE  GENERATOR. 

C (CF.  BOX  ANO  MULLER,  ANN.  MATH.  STAT.,  VOL.  28,  PP.  610-613, 
C 1958) 

C THE  INPUT  PARAMETERS  ARE 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


NFC  = NUfiee^  OF  FREQUENCY  COMPONENTS, 

UFS  = UNIT  Frequency  spacing  of  the  spectrum, 

IFCMN  = MINIMUM  OF  FREQU ENC Y-COMPONE NT  NUM3E«  FOR 
WHICH  COMPONENTS  OF  l/F  ARE  TO  BE  GENERATEO 
(MUST  Ut  2 OR  GREATER), 

IFCMX  = MAXIMUM  OF  FRECUE NCY -COMPONENT  NUMBER  FOR 

WHICH  COMPONENTS  OF  \/F  ARE  TO  BE  GENERATED  • 
(MUST  BE  NFC  OR  LESS) . 

THE  OUTPUT  PARAMETER  IS 

VF  = COMPLEX  ARRAY  OF  DIMENSION  NFC  WHERE  THE 
VOLTAGE  SPECTRUM  COMPONENTS  OF  THE  OUTPUT 
SIGNAL  ARE  TO  BE  STORED. 
this  SUBROUTINE  CALLS  THE  RANF  FUNCTION, 

UECLARATION  STATEMENTS 

DIMENSION  YF{b«.) 

DATA  TWOPI/6.2831!JdS071F939/ 

DATA  LUN/6/ 

PRELIMINARY  PROCESSING 
IJ  NC=NFC 


U0=UF3 


IMN=IFCMN 

IMX=IFCMX 


IF  (NO  .1 

LT.4) 

GO 

TO 

90 

IF (UO. 

LE.( 

).0) 

GO 

TO 

90 

IF( IMN 

.LT. 

.2) 

GO 

TO 

90 

IFdMX 

.GT, 

.NO) 

GO 

TO 

90 

IF ( IMN 

.GT, 

.IMX) 

GO 

TC 

90 

.RS  THE 

VF 

ARRAY . 

30  NT2=H0tN0 

DO  31  J=1,NT2 

YF( J)=3.U 

31  CONTINUE 

C GENERATES  A PAIR  OF  NORMAL  DEVIATES  AS  THE  REAL  AND  IMAGINARY 
C PARTS  OF  EACH  SPECTRAL  COMPONENT  OF  GAUSSIAN  NOISE  AND  STORES 
C THE  RESULTS  IN  THE  VF  ARRAY. 

AO  CF=2.J*U0 

00  41  I=IMN,IMX 
PP=RANF(0) 

A=SQRT (-CF*AL0G(PP) ) 

PH=TWOPI*RANF (0) 


J=I  + I 

VF ( J-1) =A*COS (PH) 

VF(J)  =A»SIN(PH) 

41  CONTINUE 
RETURN 
C ERROR  EXIT 

90  WRITE  (LJN,2090)  NO ,U0 , IMN , IMX 
RETURN 

C FORMAT  STATEMENTS  FOR  ERROR  MESSAGES 
2090  FORMAT  (1X-/31H  •**  IMPROPER  INPUT  VALUE(S)./ 

1 3X.5HNFC  =,I6,  5X,5HUFS  s,Elt.3, 

2 5X.7HIFCMN  = , 16 , 5X , 7HIFCMX  =,Ib/ 

3 35H  ERROR  3ETECTEO  IN  ROUTINE  SGHGN  /) 


END 
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APPENDIX  B 

USER  MANUAL  OF  THE  MODEL 


Names  and  brief  descriptions  of  all  subprograms  included  in  the  model 
are  listed  in  table  B-1,  which  is  a reproduction  of  table  1 given  in  section 
4 of  the  tex,t.  Table  B-1  is  followed  by  user  write-ups  of  all  the  subpro- 
grams in  alphabetical  order. 


Table  B-1.  List  of  Subroutines  Included  in  the  Model 

1 '■  - 1 - 1 --  1 ■ 1 

Name 

Brief  Description  of  Action 

BPFCH 

Simulates  a bandpass  filter  of  Chebyshev  type. 

CBCCT2 

Simulates  an  RF  combining  circuit  that  linearly  combines 
RF  signals. 

two 

CBCCT3 

Simulates  an  RF  combining  circuit  that  linearly  combines 
RF  signals. 

three 

DEMAM 

Simulates  an  AM  demodulator. 

DEMFM 

Simulates  an  FM  demodulator. 

DEMFSK 

Simulates  an  FSK  demodulator. 

DEMPHM 

Simulates  a phase-modulation  demodulator. 

DEMSSB 

Simulates  an  SSD  demodulator. 

OFT 

Performs  discrete  Fourier  transform. 

FRBPFC 

Calculates  frequency  response  of  a bandpass  filter  of 
Chebyshev  type. 

MODAM 

Simulates  an  AM  modulator. 

MODFM 

Simulates  an  FM  modv^lator. 

MODFSK 

Simulates  an  FSK  modulator. 

MODPHM 

Simulates  a phase-modulation  modulator. 

MODSSB 

Simulates  an  SSB  modulator. 

PSPHCT 

Simulates  a psophometer  for  commercial -telephone  circuit. 

PSPHPT 

Simulates  a psophometer  for  program-transmission  circuit. 

SGPLS 

Simulates  a pulse-signal  generator. 

SGWGH 

Simulates  a white-Gaussian-Noise  generator. 

* 

I 

1 
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BPFCH 


BPFCH 


Purpose;  This  subroutine  simulates  a bandpass  filter  of  Chebyshev  type. 


Fortran  Calling  Statement: 

CALL  BPFCH  (NFC,UFS,N0B,RPB,FCB,BWB,VF1,VF2) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  positive  but 
not  greater  than  2048), 

UFS  = unit  frequency  spacing  (must  be  positive), 

nOb  = order  of  the  bandpass  filter  (must  not  be  ne'^ative  or 
greater  than  100), 

RPB  = ripple  (dB)  within  the  passband  of  the  bandpass  filter 
(must  he  nonnegative), 

FCB  = center  frequency  of  the  bandpass  filter  (must  be 
positive), 

BWB  = 3-dB  bandwidth  of  the  bandpass  filter  (must  be 
positive) , 

VFl  = complex  array  of  dimension  NFC  containing  the  spec- 
trum components  of  the  input  voltage, 

and  the  output  parameter  is 

VF2  = complex  array  of  dimension  NFC  where  the  spectrum 
components  of  the  output  voltage  are  to  be  stored. 

Error  Message;  When  this  subroutine  is  called  with  improper  input  values,  an 
error  message  will  be  written  on  a standard  output  unit,  and  a normal 
exit  will  be  taken  from  this  subroutine. 

Note:  This  subroutine  calls  the  FRBPFC  subroutine. 
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CBCCT2 


CBCCT2 


Purpose:  This  subroutine  simulates  a radio-frequency  (RF)  combining  circuit 

that  linearly  combines  two  RF  signals  with  a specified  ratio.  The  level 
of  the  first  input  signal  is  unchanged  in  this  circuit.  • 

Fortran  Calling  Statement: 

CALL  CBCCT2  (NFC,VF1,VF2,RDB,VF3) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  in  the  spectrum  of 
each  signal, 

VFl  = complex  array  of  dimension  NFC  containing  the  ventage 
spectrum  components  of  the  first  input  signal, 

VF2  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  second  input  signal, 

RDB  = signal-to-interference  ratio  or  signal-to-noise- 
density  ratio  between  the  first  and  the  second 
input  signals  in  dB, 

and  the  output  parameter  is 

VF3  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  composite  signal  are 
to  be  stored. 


Error  Message: 


None. 


CBCCT3 


CBCCT3 


Purpose:  This  subroutine  simulates  a radio-frequency  (RF)  combining  circuit 

that  linearly  combines  three  RF  signals  with  specified  ratios.  The  level 
of  the  first  input  signal  is  unchanged  in  this  circuit.  ■ 


Fortran  Calling  Statement: 

CALL  CBCCT3  (NFC,VF1,VF2,VF3,R12DB,R13DB,VF4) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  in  the  spectrum  of  each 
signal , 

VFl  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  first  input  signal, 

VF2  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  second  input  signal, 

VF3  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  third  input  signal, 

R12DB  = signal-to-interference  ratio  or  signal -to-noise-density 
ratio  between  the  first  and  the  second  input  signals  in 
dB, 

R13DB  = signal-to-interference  ratio  or  signal -to-noise-density 
ratio  between  the  first  and  the  third  input  signals  in 
dB, 

and  the  output  parameter  is 

VF4  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  composite  signal  are 
to  be  stored. 

Error  Message:  None. 
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Fortran  Calling  Statement: 

CALL  DEMAM  (NFC.VF.SF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2,3 14), 

VF  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  demodulator  input  signal, 

and  the  output  parameter  is 

SF  = complex,  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  demodulated  signal  are 
to  be  stored. 


Error  Message:  None. 


Note;  This  subroutine  calls  the  DFT  subroutine. 
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DEMFM 


DEMFM 


Purpose:  This  subroutine  simulates  a frequency-modulation  (FM)  demodulator. 

Fortran  Calling  Statement: 

CALL  DEMFM  (NFC,UFS,VF.FCD,SF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2,3,. ...14), 

UFS  = unit  frequency  spacing  of  the  spectrum, 

VF  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  demodulator  input  signal, 

FCD  = center  frequency  of  the  demodulator, 
and  the  output  parameter  is 

SF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  demodulated  signal  are 
to  be  stored. 

Error  Message:  None. 

Note:  This  subroutine  calls  the  DFT  subroutine. 


63 
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DEMFSK 


DEMFSK 


Purpose:  This  subroutine  simulates  a frequency-shift-keying  (FSK)  demodulator. 


Fortran  Calling  Statement: 

CALL  DEMFSK  (NFC,UFS,VF,FCD,NIB,IB) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2,3,. ...11), 

UFS  = unit  frequency  spacing  of  the  spectrum, 

VF  = complex  array  of  dimension  MFC  containing  the  voltage 
spectrum  components  of  the  demodulator  input  signal, 

FCD  = center  frequency  of  the  demodulator, 

NIB  = number  of  information  bits  in  a period, 
and  the  output  parameter  is 

IB  = array  of  dimension  NIB  where  the  demodulated  informa- 
tion bits  are  to  bo  stored. 

Error  Message:  When  NFC  exceed  2048  (=2**11),  an  error  message  will  be  writ- 
ten on  a standard  output  unit,  and  a normal  exit  will  be  taken  from  this 
subroutine. 


Note:  This  subroutine  calls  the  DFT  subroutine. 
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DEMPHM 


DEMPHM 


Purpose:  This  subroutine  simulates  a phase-modulation  demodulator. 


Fortran  Calling  Statement: 

CALL  DEMPHM  (NFC.UFS.VF.SF) 

where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 
M = 2,3 14), 

UFS  = unit  frequency  spacing  of  the  spectrum. 


VF  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  demodulator  input  signal. 


and  the  output  parameter  is 


SF  * complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  demodulated  signal  are 
to  be  stored. 


The  output,  SF,  is  normalized  in  such  e way  that  its  power  will 
equal  unity  when  the  input  signal  is  phase  modulated  by  a single 
sinusoid  with  a modulation  index  of  unity. 


Error  Message:  None. 


Note:  This  subroutine  calls  the  DFT  subroutine. 


DEMSSB 


DEMSSB 


Purpose:  This  subroutine  simulates  a single-sideband  (SSB)  demodulator. 

Fortran  Calling  Statement: 

CALL  DEMSSB  (NFC,UFS,VF,FRD,SF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2. 3,..., 14), 

UFS  = unit  frequency  spacing  of  the  spectrum, 

VF  = complex  array  of  dimension  NFC  containing  the  voltage 
soectrum  components  of  the  demodulator  input  signal, 

FRD  = reference  frequency  of  the  demodulator, 
and  the  output  parameter  is 

SF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  demodulated  signal  are 
to  be  stored. 

Error  Message:  None. 

Note:  This  subroutine  calls  the  RANF  function. 
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DFT 


DFT 


Purpose:  This  subroutine  performs  discrete  Fourier  transform  between  time 

and  frequency  domains. 

Fortran  Calling  Statement: 

CALL  DFT  (N,X1,ISGN,X2) 

where  the  input  parameters  are 

M = dimension  of  the  XI  and  X2  arrays, 

XI  = complex  array  of  dimension  N containing  the  input 
sequence  to  be  transformed, 

ISGN  = sign  in  the  exponent 

= +1  for  transform  from  frequency  domain  to  time  domain 
= -1  for  transform  from  time  domain  to  frequency  domain, 
and  the  output  parameter  is 

X2  = complex  array  of  dimension  N,  where  the  transformed 
sequence  is  to  be  stored. 

Restriction;  N must  be  equal  to  2**M,  where  M is  an  integer  from  0 to  14 
inclusive. 

Error  Message:  If  N violates  the  above  restriction,  an  error  message  will  be 

written  on  the  standard  output  unit,  and  a normal  exit  will  be  taken 
from  the  subroutine. 

Timing:  Approximately  0.1  N milliseconds  (on  the  CDC-6600  computer). 

Note:  This  is  a modified  version  of  the  discrete  complex  fast  Fourier  trans- 

form package  programmed  by  L.  David  Lewis  and  Marie  West  of  the  Space 
Environment  Laboratory,  National  Oceanic  and  Atmospheric  Administration, 
U.S.  Department  of  Commerce,  Boulder,  Colorado.  Their  subroutine  pack- 
age is  a Fortran  translation  from  an  Algol  procedure  written  by  R.  C. 
Singleton,  "Algorithm  339,  An  Algol  procedure  for  the  fast  Fourier  trans- 
form with  arbitrary  factors,"  Common.  ACM,  vol.  11,  no.  11  (Nov.  1968), 
pp.  776-779,  and  vol.  12,  no.  3 (March  1969),  p.  189. 
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FRBPFC 


FRBPFC 


Purpose:  This  subroutine  calculates  frequency  response  of  a Chebyshev-type 
bandpass  filter,  (cf.  J.  E.  Storer,  Passive  Network  Synthesis,  McGraw- 
Hill,  ch.  30) 

Fortran  Calling  Statement: 

CALL  FRBPFC  (NO,RPL,FC,BW,UFS,IFCtIN,IFCMX,HF) 

where  the  input  parameters  are 

NO  = order  of  the  bandpass  filter  (must  be  nonnegative 
but  not  greater  than  100), 

RPL  = ripple  (dB)  within  the  passband  (must  be  nonnegative), 

FC  = center  frequency  (must  be  positive), 

BW  = 3-dB  bandwidth  (must  be  positive), 

UFS  = unit  frequency  spacing  (must  be  positive), 

IFCMN  = minimum  of  frequency-component  number  for  which 
frequency  response  is  to  be  calculated, 

IFCMX  = maximum  of  frequency-component  number  for  which 
frequency  response  is  to  be  calculated, 

and  the  output  parameter  is 

HF  = complex  array  of  dimension  IFCMX  or  greater  where  the 
calculated  frequency  response  for  0FS*(I-1)  is  to  be 
stored  as  its  Ith  element. 

Error  Message:  When  this  subroutine  is  called  with  improper  input  values,  an 

error  message  will  be  written  on  a standard  output  unit,  and  a normal 
exit  will  be  taken  from  this  subroutine. 
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MODAM 


MODAM 


Purpose:  This  subroutine  simulates  an  amplitude-modulation  (AM)  modulator. 


Fortran  Calling  Statement: 

CALL  MODAM  (NFC,UFS,SF,FC,FLMI,VF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2.3.. ..,14), 

UFS  = unit  frequency  spacing  of  the  spectrum  (must  be 
positive), 

SF  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  modulating  signal  (must 
include  at  least  one  nonzero  element), 

FC  = center  frequency  of  the  modulated  signal  (must  be 
smaller  than  (NFC-1)*UFS) , 

FLMI  = modulation  index  (must  not  exceed  1.0), 
and  the  output  parameter  is 

VF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  modulated  signal  are 
to  be  stored. 

Error  Message:  When  this  subroutine  is  called  with  improper  input  values,  an 

error  message  will  be  written  on  a standard  output  unit,  and  a normal 
exit  will  be  taken  from  this  subroutine. 


Note:  This  subroutine  calls  the  RANF  function. 
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MODFM 


MODFM 


Purpose:  This  subroutine  simulates  a frequency-modulation  (FM)  modulator; 

Fortran  Calling  Statement: 

CALL  MODFM  (NFC,UFS,SF,FC,FDV,VF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2,3,. ..,14), 

UFS  = unit  frequency  spacing  of  the  spectrum  (must 
be  positive), 

SF  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  modulating  signal  (must 
include  at  least  one  nonzero  element), 

FC  = center  frequency  of  the  modulated  signal  (must  be 
smaller  than  (NFC-1)*UFS) , 

FDV  = frequency  deviation, 
and  the  output  parameter  is 

VF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  modulated  signal  are 
to  be  stored. 

Error  Message:  When  this  subroutine  is  called  with  improper  input  values,  an 

error  message  will  be  written  on  a standard  output  unit,  and  a normal 
exit  will  be  taken  from  this  subroutine. 


Note:  This  subroutine  calls  the  DFT  subroutine  and  the  RAMF  function. 
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MODFSK  NODFSK 

Purpose:  This  subroutine  simulates  a frequency-shift-keying  (FSK)  modulator. 


Fortran- Calling  Statement: 

CALL  MODFSK  (NFC,UFS,NIB,IB,FC,FS,VF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2,3,. ..,14), 

UFS  = unit  frequency  spacing  of  the  spectrum  (must  be 
positive) , 

NIB  = number  of  information  bits  in  a period  (must  not 
exceed  NFC/2), 

IB  = integer  array  of  dimension  NIB  containing  the 
information  bits, 

FC  = center  frequency  of  the  modulated  signal  (must  be 
smaller  than  (NFC-1)*UFS) , 

FS  = frequency  shift  (must  be  smaller  than  2*FC  and 
2*(NFC*UFS-FC)), 

and  the  output  parameter  is 

VF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  modulated  signal  are 
to  be  stored. 

Error  Message:  When  this  subroutine  is  called  with  improper  input  values,  an 

error  message  will  be  written  on  a standard  output  unit,  and  a normal 
exit  will  be  taken  from  this  subroutine. 

Note:  This  subroutine  calls  the  DFT  subroutine  and  the  RANF  function. 
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MODPHM 


MODPHH 


Purpose:  This  subroutine  simulates  a phase-modulation  modulator. 


Fortran  Calling  Statement: 

CALL  MODPHM  (NFC.UFS ,SF,FC,FLMI ,VF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

K = 2,3,, ..,14), 

UFS  = unit  frequency  spacing  of  the  spectrum  (must 
be  positive), 

SF  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  modulating  signal  (must 
include  at  least  one  nonzero  element), 

FC  = center  frequency  of  the  modulated  signal  (must  be 
smaller  than  (NFC-1)*UFS) , 

FLMI  = modulation  index  in  radians, 
and  the  output  parameter  is 

VF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  modulated  signal  are 
to  be  stored. 


Error  Message:  When  this  subroutine  is  called  with  improper  input  values, 

an  error  message  will  be  written  on  a standard  output  unit,  and  a 
I normal  exit  will  be  taken  from  this  subroutine. 

t 

}■ 

i Mote:  This  subroutine  calls  the  DFT  subroutine  and  the  RANF  function. 


MODSSB 


MODSSB 


Purpose:  This  subroutine  simulates  a single-sideband  (SSB)  modulator. 


Fortran  Calling  Statement: 

CALL  MODSSB  (NFC,UFS,SF,FR,VF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2,3.. ..,14), 

UFS  = unit  frequency  spacing  of  the  spectrum  (must 
be  positive), 

SF  = complex  array  of  dimension  NFC  containing  the  voltage 
spectrum  components  of  the  modulating  signal  (must 
include  at  least  one  nonzero  element), 

FR  = reference  frequency  of  the  modulated  signal  (must  be 
smaller  than  (NFC-1)*UFS) , 

and  the  output  parameter  is 

VF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  modulated  signal  are 
to  be  stored. 

Error  Message:  When  this  subroutine  is  called  with  improper  input  values, 

an  error  message  will  be  written  on  a standard  output  unit,  and  a nor- 
mal exit  will  be  taken  from  this  subroutine. 


Note:  This  subroutine  calls  the  RANF  function. 


PSPHCT 


PSPHCT 


► Purpose:  This  subroutine  simulates  a psophometer  for  commercial -telephone 

i circuits,  (cf.  CCITT  Recommendation,  p.  53,  1973) 

t 

I Fortran  Calling  Statement: 

CALL  PSPHCT  (NFC,UFS,SF,PW) 

' where  the  input  parameters  are 

. NFC  = number  of  frequency  components, 

UFS  = unit  frequency  spacing  of  the  spectrum  in  hertz, 

SF  = complex  array  of  dimension  MFC  containing  the 
i voltage  spectrum  components  of  the  input  signal, 

and  the  output  parameter  is 

PW  = psophometrically-weighted  signal  power. 


Error  Message:  None. 


i 
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PSPHPT 


PSPHPT 


Purpose:  This  subroutine  simulates  a psophometer  for  program-transmission 

circuits,  (cf.  CCITT  Recommendation,  p.  [13,  1973) 

Fortran  Calling  Statement: 

CALL  PSPHPT  (NFC,UFS,SF,PW) 

where  the  input  parameters  are 

NFC  = number  of  frequency  components, 

UFS  = unit  frequency  spacing  of  the  spectrum  in  hertz, 

SF  = complex  array  of  dimension  NFC  containing  the 

voltage  spectrum  components  of  the  input  signal, 

and  the  output  parameter  is 

PW  = psophometrically-weighted  signal  power. 

Error  Message:  None. 
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Purpose:  This  subroutine  simulates  a pulse-signal  generator. 


Fortran  Calling  Statement: 

CALL  SGPLS  (NFC,UFS,FC,PRF,TR,TW,TF,DVF,VF) 
wherf  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

M = 2,3,. ..,14), 

UFS  = unit  frequency  spacing  of  the  spectrum  (must  be 
positive), 

FC  = center  frequency  of  output  signal  (must  be 
smaller  than  NFC*UFS), 

PRF  = pulse  repetition  frequency, 

TR  = pulse  rise  time, 

TW  = pulse  width  at  the  top, 

TF  = pulse  fall  time, 

DVF  = maximum  frequency  deviation  (must  be  smaller 
than  FC  and  NFC*UFS-FC), 

and  the  output  parameter  is 

VF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  output  signal  are  to 
be  stored. 

Error  Message:  When  this  subroutine  is  called  with  improper  input  values, 

an  error  message  will  be  written  on  a standard  output  unit,  and  a nor- 
mal exit  will  be  taken  from  this  subroutine. 

Note:  This  subroutine  calls  the  DFT  subroutine  and  the  RANF  function. 
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SGWGN 


SGWGN 


Purpose:  This  subroutine  simulates  a white-Gaussian-noise  generator. 


Fortran  Calling  Statement: 

CALL  SGWGN  (NFC,UFS,IFCMN,IFCMX,VF) 
where  the  input  parameters  are 

NFC  = number  of  frequency  components  (must  be  2**M, 

* M = 2,3,. ...14), 

UFS  = unit  frequency  spacing  of  the  spectrum  (must  be 
positive) , 

IFCMN  = minimum  of  frequency-component  number  for  which 
components  of  VF  are  to  be  generated  (must  be 
2 or  greater), 

IFCMX  = maximum  of  frequency-component  number  for  which 
components  of  VF  are  to  be  generated  (must  be 
IFCMN  or  greater,  but  not  greater  than  NFC), 

and  the  output  parameter  is 

VF  = complex  array  of  dimension  NFC  where  the  voltage 
spectrum  components  of  the  output  signal  are  to 
be  stored. 

Error  Message:  When  this  subroutine  is  called  with  improper  input  values,  an 

error  message  will  be  written  on  a standard  output  unit,  and  a normal 
exit  will  be  taken  from  this  subroutine. 


Note:  This  subroutine  calls  the  RANF  function. 
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